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ABSTRACT 



F. Muller Sanchez 1 , A. Sternberg 4 



In a sample of local active galactic nuclei (AGNs) studied at a spatial resolution on the order of 
10 pc, we show that the interstellar medium traced by the molecular hydrogen v—1-0 S(l) line at 2.1 
/im forms a geometrically thick, clumpy disk. The kinematics of the molecular gas reveals general 
rotation, although an additional significant component of random bulk motion is required by the high 
local velocity dispersion. The size scale of the typical gas disk is found to have a radius of ^30 pc with 
a comparable vertical height. Within this radius the average gas mass is estimated to be ~10 7 M 
based on a typical gas mass fraction of 10%, which suggests column densities of Nh ~ 5 x 10 23 cm~ 2 . 
Extinction of the stellar continuum within this same region suggests lower column densities of Nh ~ 
2 x 10 22 cm -2 , indicating that the gas distribution on these scales is dominated by dense clumps. 
In half of the observed Seyfert galaxies this lower column density is still great enough to obscure the 
AGN at optical/infrared wavelengths. We conclude, based on the spatial distribution, kinematics, and 
column densities that the molecular gas observed is spatially mixed with the nuclear stellar population 
and is likely to be associated with the outer extent of any smaller scale nuclear obscuring structure. 
Furthermore, we find that the velocity dispersion of the molecular gas is correlated with the star 
formation rate per unit area, suggesting a link between the two phenomena, and that the gas surface 
density follows known "Schmidt-Kennicutt" relations. The molecular/dusty structure on these scales 
may be dynamic since it is possible that the velocity dispersion of the gas, and hence the vertical disk 
height, is maintained by a short, massive inflow of material into the nuclear region and/or by intense, 
short-lived nuclear star formation. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: 
Seyfert — infrared: galaxies 
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1. INTRODUCTION 

Although there is a great deal of evidence supporting 
the widely accepted active galactic nuclei (AGNs) uni- 
ficati on model (e.g. lAntonu cci 1993, lUrry fc PadovanH 
Il995f ) , characterization, or even proof, of the key model 
component, an obscuring medium, has been difficult to 
attain. A popular version of the unification model pro- 
poses the existence of an optically and geometrically 
thick component of gas/dust that obscures the optical 
and ultraviolet (UV) AGN emission along some fraction 
of the lines of sight such that, depending on the view- 
ing angle, the observed characteristics differ, resulting 
in type 1 and 2 AGNs. In this scenario, the obscur- 
ing medium, which is classically thought of as a static 
canonical "torus" surrounding the AGN accretion disk, 
is thought to begin at the dust sublimation radius and 
extend outward. 

* Based on observations at the ESO Very Large Telescope (60. A- 
9235, 070.B-0649, 070.B-0664, 074.B-9012, 076.B-0098, 076.B- 
0117, 077.B-0514). 

' Based on observations at the W. M. Keck Observatory, which is 
operated as a scientific partnership among the California Institute 
of Technology, the University of California and the National Aero- 
nautics and Space Administration. The Observatory was made 
possible by the generous financial support of the W.M. Keck Foun- 
dation. 

1 Max Planck Institut fur extraterrestrische Physik, Postfach 
1312, 85741, Garching, Germany 

2 Department of Physics and Astronomy, University of Califor- 
nia, Los Angeles, CA, 90095-1562, United States 

3 Department of Physics, 366 Le Conte Hall, University of Cali- 
fornia, Berkeley, CA, 94720-7300, United States 

4 School of Physics and Astronomy, Tel Aviv University, Tel Aviv 
69978, Israel 



Thermal dust emission has been spatially resolved 
around AGNs by infrare d (IR) interferometric techniques 
in the AGN NGC 1068 tiJaffe et alJl200l IWeieelt et all 
12004 IPoncelet et all 120061 IPoncelet et all l2007t) and 
Circinus (|Tristram et al.ll2007f ). This emission is inter- 
preted as reradiated emission from dust heated by the 
optical/ultraviolet emission originating in the accretion 
disk. The data suggest a two-component structure in 
both galaxies consisting of a subparsec, warm, disklike 
component that is surrounded by a slightly lower tem- 
perature, geometrically thick, torus-like component that 
extends out to a radius of a few parsecs. In several 
nearby AGNs H2 O masers are found to trace a subparsec 
warped disk (e.g. iHerrnstein et al.lll996l iGreenhill et al.1 
[l996llGreenhill et al.ll2003[ ). and in NGC 1068 and Circi- 
nus these warped disks have spatial scales and orienta- 
tions that are consistent with the warm disklike com- 
ponent detected in the IR observations. The compact 
size of the thermal dust emission implies that this ob- 
scuring medium is in a clumpy distribution, rather than 
smo othly d i stribut ed as assumed in early torus mod- 
els (jElitzurl l200 6b). In obscured nuclei, X-ray obser- 
vations imply column densities on these parsec scales 
of Nh ~ 10 22 cm -2 to greater than Nh ~ 10 24 cm -2 
(|Risaliti et all I1999L ITreister et aLI 120041 ). which is in- 
deed high enough to obscure the optical and ultraviolet 
emission from the accretion d isk (>10 22 cm~ 2 for Seyfert 
galaxies; ITreister et alll2004 . 

Models of clumpy tori predict that the obscuring 
region extends from subparsec radii out to sizes on 
the orde r of tens of parsecs i n typical Seyfert galax- 
ies (e.g. ICameron et al.l Il9 93. Schartm ann et al. | 120051 . 
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TABLE 1 

Summary of AGN data included in this study 



IFU AO Observations 



ISAAC Observations 



Object 


Classification a 


D(Mpc) a 


pc/" 


Band 


Res. b (") 


Rcs. b (pc) 


Date 


PAs (°) 


Date 










SINFONI 












NGC 1097 


LINER, Sy 1 


18 


88 


H+K 


0.25 


22.0 


2005 Oct 


-33,57 


2006 Sep 


NGC 3227 c 


Syl 


17 


80 


K 


0.09 


7.2 


2004 Dec 


-30,60 


2006 Jan 


NGC 3783 


Syl 


12 


200 


K 


0.18 


36.0 


2005 Mar 






NGC 7469 


Syl 


66 


320 


K 


0.15 


48.0 


2004 Jul 






Circinus c 


Sy2 


4 


20 


K 


0.22 


4.4 


2004 Jul 


-65,25 


2006 Jan,Feb 


NGC 1068 c 


Sy2 


14 


70 


H+K 


0.08 


5.6 


2005 Oct, 2006 Nov 














OSIRIS 












NGC 3227 


Syl 


17 


80 


K 


0.07 


5.6 


2006 Mar, 2008 Jan 






NGC 4051 


Syl 


10 


18 


K 


0.12 


5.8 


2006 Apr, 2008 Jan 






NGC 4151 


Syl 


11 


68 


K 


0.11 


7.5 


2006 Mar 






NGC 6814 


Syl 


22 


107 


K 


0.17 


18.2 


2006 Apr, Sep 






NGC 7469 


Syl 


67 


320 


K 


0.11 


35.2 


2006 Sep 







a Classification and distance (D) arc from D07 and the NASA/IPAC Extragalactic Database. 

b Spatial resolution (FWHM) estimated from the data itself, using the methods described in § 2.2. 

c References to detailed s tudies of individual objects: Circinus JMiiller Sanchez et ail 120061) . NGC 3227 JDavies et al.M2006D . NGC 1068 
IMiiller Sanchez c t al.ll2009l) . 



Honig et alj 120061 iFritz et all 120061 ISchartmann et all 
20081 Schartmann et al.ll2009f ) . Little is known about the 
characteristics of the gas and dust on these scales beyond 
the relatively small-scale structures inferred from the IR 
interferometry and H2O masers measurements. Obser- 
vations of the thermal dust emission are temperature 
dependent and the inferred radii are thus temperature- 
weighted. These measurements therefore represent only 
the warm dust and do not trace any potentially cooler 
extended structures in the nuclear region. Due to the 
glare from the AGN, studies probing the obscuring re- 
gion predicted by models are challenging, and those 
that have been carried out focus primarily on Seyfert 
2 ga laxies to reduce the impact of the AGN e mission 
(e.g. ICid Fernandes et al1l2004 ISilge et alj|2005| ). How- 
ever, with adaptive optics (AO) the AGN emission can 
now be confined to a smaller region, making possible a 
detailed study of the nuclear region (central — 100 pc) 
in even bright Seyfert 1 galaxies. This provides a more 
direct view of the nuclear region and eliminates poten- 
tial biases imposed by the greater amount of obscuration 
present in Seyfert 2s. It is therefore now possible to es- 
tablish the respective roles of the smaller and larger scale 
components in fulfilling the functions, such as obscura- 
tion of the AGN, of a unifying torus. 

In an effort to better characterize the obscuring 
medium present out to the several tens of parsec scales 
predicted by models, and investigate its association with 
the nuclear starburst activity and fueling of the AGN, 
we have observed a sample of nine local AGNs using K- 
band AO assisted integral field spectroscopy. The bright 
molecular hydrogen emission from the 1-0 S(l) rovibra- 
tional transition at 2.1218 /zm is used to measure, on 
scales of —10 pc, the two-dimensional (2-D) distribu- 
tion and kinematics of the nu clear gas within the cen - 
tral -100 pc. As reported by IChristopher et alj (|2005t ) 
to be the case in the Galactic Center, this warm molec- 
ular gas traces the cold/dense gas within the nuclear re- 
gion. Sample selection, observations, and data reduction 
are discussed in § [2] and the general properties of the 
molecular gas are described in § [3l Presented in § |4] 



are correlations between the observed properties of the 
molecular gas and the inferred nuclear star formation. 
Possible scenarios for maintaining the estimated vertical 
structure of the gas disk are investigated in § [5j and in 
§ [6] the role of the molecular hydrogen on these scales in 
obscuring the AGN and its connection with the nuclear 
starburst activity are discussed. The overall conclusions 
of the study are outlined in § \7\ while a more detailed 
discussion of each of the observed galaxies can be found 
in the Appendix. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Sample Selection 

The nine galaxies were observed at the ESO Very Large 
Telescope (VLT) and the W. M. Keck Observatory. Al- 
though the AGNs observed were not selected to comprise 
a complete sample, both Seyfert Is and 2s, as well as a 
LINER (low-ionization nuclear emission-line region), are 
included. Table[T]lists the full sample of galaxies selected. 

The criteria for target selection of the galaxies observed 
at the VLT were (1) the nucleus is bright enough to use 
for AO correction, (2) the galaxy is close enough that 
small spatial scales can be resolved (—20 pc), and (3) the 
galaxies are well studied such that complementary data 
can be found in the literature (e.g. radio, submillimeter 
observations). Although, the first criterion biases the 
sample to type 1 Seyferts, two of the six galaxies observed 
are Seyfert 2s. 

The sample observed at the Keck Observatory was se- 
lected from all AGNs with a black hole mass estimate 
from reverberation mapping and a redshift such that the 
2.1218 /jm molecular hydrogen emission and stellar 2.3 
/xm CO bandheads are observable in the K band. This 
sample is therefore composed of relatively nearby Seyfert 
1 galaxies, all of which have a nucleus that is bright 
enough for either a correction with the natural guide 
star (NGS) AO system or as a tip-tilt guide to be used 
in conjunction with the laser guide star (LGS) AO sys- 
tem. Both the OSIRIS and SINFONI subsamples contain 
NGC 3227 and NGC 7469, and a comparison of the data 
from each instrument is presented in the Appendix. 
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2.2. Integral Field Spectra from SINFONI 

Integral field spectroscopy was obtained for several 
AGN at the VLT UT4 with the AO near-IR inte- 
gral field spectrogra ph SINFONI (jEisenhauer et al.ll2003l . 
iBonnet et all 12004 ) . The pixel scale of the R ~ 4300 
if-band spectra (approximately 1.95-2.45 /zm) is either 
0". 125x0". 25 or 0".05x0".l, which results in a field of 
view (FOV) of 0".8x0".8 or 3.2"x3.2", respectively. For 
each galaxy, the nucleus was used as the reference for 
a near-diffraction-lim ited correction by the AO module 
(|Bonnet et al.1 l2003h . The exposure times for Circinus, 
NGC 3783, and NGC 1068 were 30, 40, and 110 minutes, 
respectively, and for the remaining galaxies 60 minutes of 
on-source data were obtained. A summary of the data, 
including the specific spatial resolution achieved for each 
galaxy, is given in Table [1] 

The data were reduced using the SINFONI custom re- 
duction package SPRED ([Abuter et al.ll2006f ). which in- 
cludes all of the typical reduction steps applied to near- 
IR spectra with the additional routines necessary to re- 
construct the data cube. After background subtraction, 
which is performed with sky frames obtained interspersed 
with the on-source exposures, the data are flat-fielded 
and corrected for dead/hot pixels. Telluric correction 
and flux calibration are carried out using A- and B-type 
stars. Residuals from the OH lin e emiss i on are minimized 
using the methods outlined in iDaviesI (|2007l ). The flux 
calibration is confirmed by comparison to 3" aperture 
Two Micron All Sky Survey and l"-3" aperture NACO 
or Hubble Space Telescope {HST) NICMOS images and 
was found to be consistent to within 20%. 

The spatial resolution of the SINFONI data (Table 
[1]) is determined from the data themselves rather than 
from additional point-spread function (PSF) calibration 
frames. This is done because of the uncertainty involved 
in using a standard star to recreate the AO performance 
achieved with an AGN which has spatially extended 
background galaxy light, as well as to avoid mismatches 
due to the varia bility of the ambient seeing. Two meth- 
ods, outlined in iDavies et all l2004al are used to deter- 
mine the PSF directly from the data, which also has the 
advantage of including all effects associated with the con- 
struction of the data cube. The first method makes use 
of the broad Br7 (A2.1661 /im) emission assumed to be 
coming from the broad- line region (BLR), which is ex- 
pected to be unresolv ed for the galaxies in this sample 
(jPeterson et all 12004 ). The second method applies the 
fact that the intrinsic equivalent width of the CO 2-0 
2.29 /x m bandhead is expected to fall within a limited 
range (jOliva et al.lll995l iForster Schreiberll2000L Davies 
et al. 2007, hereafter D07), and therefore the observed 
Wco2-o can be used to determine the dilution of the stel- 
lar feature due to the AGN continuum emission, which, 
at this resolution is also a point source (predicted size s 
are less than 2 pc; Uaffe et al.ll2004 iTiTstram et al.|[2007D . 

An analysis of the nuclear star formation properties 
(star formation rates (SFRs), time since the last star- 
forming episode, stellar distribution, etc.) in each of 
the galaxies in this subsample is presented in D07. 
Detailed studies of a few of the galaxies in the sub- 
sample can also be found in the literatu re: Circinus 
(iMuller Sanchez et all I2006T). NGC 3227 (IDavies et all 
2006j), and NGC 1068 (|Miiller Sanchez et al.ll200'9F ~ 



2.3. Integral Field Spectra from OSIRIS 

At the W. M. Keck Observatory, the in tegral field 
spectrometer OSIRIS (|Larkin et all |2006() on Keck 
II was used with both the NGS and LG S AO 
systems (Ivan Dam et all 12004 Ivan Dam et all [2006, 
IWizinowich et al.ll2006l ) Jf-band R ~ 3000 spectra, cov- 
ering 1.965-2.381 /im, were obtained with a pixel scale of 
0".035 pixel" 1 over an FOV of 0".56 x 2". 24 for each of 
the galaxies. Typically 10 minute exposures were taken 
with off-source sky exposures obtained interspersed with 
the on-source frames. Total on-source exposures vary 
from 40 minutes for NGC 4051, to 60 minutes for NGC 
4151, and 70 minutes for the remaining three galaxies. 
Of the 70 minutes of on-source exposure for NGC 7469, 
20 minutes of this was obtained with the LGS AO sys- 
tem, while the rest of the data for this galaxy, and the 
others, was obtained with a correction from the NGS 
AO system using the AGN nucleus as the reference. The 
spatial resolution achieved for each galaxy, via the meth- 
ods discussed in the previous section, are given in Table 
[TJ As a consistency check, the PSF full-width at half- 
maxima (FWHMs) estimated using these methods were 
compared to that of the PSFs determined from obser- 
vations of standard stars (using the same AO setting as 
used for the galaxies), and they are found to be con- 
sistent to within the uncertainties discussed in § 12.21 of 
standard star PSF measurements (e.g. changes in seeing 
and difficulty in reproducing the AO performance of the 
galaxy observations). 

Reduction of the OSIRIS data was carried out with the 
OSIRIS final data reduction pipeline (DRP), which per- 
forms reduction steps typically applied to near-IR spec- 
tra as well as reconstruction of the data cube. These 
steps include background subtraction, flat fielding, and 
correction for dead/hot pixels, as well as correction of 
detector nonlinearity, removal of detector cross talk, and 
wavelength calibration. Telluric correction was done us- 
ing A-type stars and was applied to the data in the DRP. 
Flux calibration of the OSIRIS data is expected to be ac- 
curate to about 30%. 

2.4. Wide FOV Long Single Slit Spectra from ISAAC 

To assess the kinematic properties of the galaxies on 
a larger scale than is measured here, long-slit data were 
also obtained with ISAAC on the VLT UT1 for three 
of the galaxies measured with SINFONI. Seeing-limited 
data were obtained as part of a service observing pro- 
gram carried out in 2006. The short wavelength arm 
of ISAAC was used in the medium resolution mode (R 
~ 3000) with a 0.122 /im narrowband filter centered at 
2.144 /zm to obtain spectra of the H2 1-0 S(l) emission 
line. Two orthogonal slit positions were obtained for each 
of the galaxies (sec Table [T] for details) with a 2' long and 
0".3 wide slit. The orientations of the two slit positions 
for each object were selected to align with the major and 
minor axes as determined from 2.2 /im HST NICMOS 
images. Exposure times totaled 40 minutes per slit posi- 
tion, and the target was nodded along the slit to perform 
sky subtraction. The ISAAC data were reduced using 
standard near-IR reduction techniques, and were carried 
out with the ESO pipeline software. 
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3. GENERAL PROPERTIES OF THE MOLECULAR GAS 

In order to draw general conclusions about the prop- 
erties of molecular hydrogen in AGNs, all of the galaxies 
in the sample are evaluated using a consistent method to 
minimize differences due to spatial resolution. The dis- 
tribution and kinematics of the molecular gas are mea- 
sured by fitting the spectrum of the H2 1-0 S(l) emission 
line with an unresolved line profile (e.g. a sky line) con- 
volved with a Gaussian, as well as a linear function to the 
line free continuum, at each spaxel. This fitting process 
gives the intrinsic velocity dispersion (i.e. instrumental 
broadening has been accounted for), and it is this value 
that is used in the subsequent analysis, figures, and ta- 
bles. The uncertainty of the fit is estimated using Monte 
Carlo techniques by refitting the best-fit Gaussian with 
added noise of the same statistics as the data. This is 
done 100 times, and the standard deviations of the best- 
fit parameters are used as the uncertainties. Typically 
the uncertainties in velocity and velocity dispersion (cr) 
are on the order of 5-15 km s _1 . The 2-D maps of the H2 
flux distribution, velocity, and a for each of the observed 
galaxies are shown in Figures [1] and O 

The gas kinematics is characterized using the method 
of kinemetry, which is a generalization of surface pho- 
tometry t o the higher order mom ents of the velocity dis- 
tribution (|Krainovic et al.l l2006). The kinematic center 
is assumed to be coincident with the nucleus defined as 
the peak in the nonstellar emission. As justified in § 13.21 
the kinematics are ultimately fit holding the position an- 
gle (PA) and inclination angle constant since neither is 
found to vary significantly in the region measured. Fur- 
thermore, the best fit PA and inclination angle of the H2 
kinematics is consistent with that fit of the stellar kine- 
matics and we conservatively assume the best-fit values 
to the stellar kinematics when possible. The azimuthal 
average of the flux distribution and kinematics in each 
galaxy are also computed, correcting for these best-fit PA 
and inclination angles, to further investigate the general 
properties of H2 in the observed AGN (Figures [3] and |4|). 

The flux distribution, rotational velocity, and velocity 
dispersion as a function of radius were also determined 
for each of the single slit ISAAC spectra. This is achieved 
by fitting a Gaussian to the emission line profile in a 
spatial region determined by binning pixels until at least 
a 5cr detection of the H2 emission is reached. A single 
Gaussian is found to be a suitable representation of the 
emission line profile in all three of the galaxies observed. 
The H2 velocity and a of the gas on the larger scales 
measured with ISAAC, as well as a comparison to the 
smaller scale kinematics measured with SINFONI, are 
shown in Fig. [5J 

We have chosen to exclude NGC 1068 from the AGN 
sample because it is the only galaxy in the sample for 
which the H2 kinematics are dominated by noncircular 
motions rather than, as for all other galaxies in the sam- 
ple, by rotation in the gravitational potential implied 
by the nuclear stellar kinematics (see § I3.2|) . The H2 
kinematics of NGC 1068 differ from circular rotation by 
as much as 100 km s" 1 , an d, as discussed in detail in 
iMMer Sanchez et all (2009), suggest an inflow of mate- 
rial toward the central BH. The high spatial resolution 
of the NGC 1068 observations can be ruled out as the 
reason for the complexity of the kinematics since the 6 



pc resolution is comparable to that achieved with four 
other galaxies (FHWM < 10 pc) and two of these have 
data with a resolution even better than that of NGC 
1068. Although we do not include NGC 1068 in our dis- 
cussion of the general H2 properties in AGNs, it is still 
presented along side the rest of the sample in all figures 
and tables for completeness. For many of the H2 proper- 
ties discussed throughout this paper, the properties mea- 
sured and derived for NGC 1068 are broadly consistent 
with those found for the rest of the sample. However, 
with the advantage of 2-D spatial data, the kinematics 
clearly reveals that the nuclear properties of NGC 1068 
are unique with respect to the other galaxies observed. 

3.1. Spatial Distribution 

As can be seen in the 2-D maps of the H2 flux distri- 
bution in Figures [1] and [2] there is a great deal of diver- 
sity in the nuclear H2 flux distributions in the sample 
galaxies. Although most of the galaxies have a relatively 
symmetric flux distribution, the 2-D maps reveal several 
nonsymmetric features, particularly in NGC 3227 and 
NGC 4151. This should not come as a surprise given 
that H2 1-0 S(l) is emitted from warm molecular gas 
(typically 1000-2000 K), and therefore the local environ- 
ment will greatly influence the observed luminosity dis- 
tribution. For this reason it should be kept in mind that 
the flux distribution is not necessarily representative of 
the actual H2 mass distribution. In fact, particularly in 
the cases of NGC 3227 and NGC 4151, in addition to 
the brighter patches of emission, H 2 is also detected at 
a significant level throughout the measured FOV. The 
H 2 flux distribution in the individual galaxies observed 
is discussed in more detail in the Appendix. 

Despite the complexity of the observed flux distribu- 
tions, the azimuthal average of the H2 emission is con- 
sistent with a half-width at half-maximum (HWHM) less 
than 35 pc in all of the AGN, with a mean HWHM for the 
sample of ~30 pc (Fig. [6] and Table [2]). Of note is that 
D07 report a nuclear stellar component with similar size 
scales in the SINFONI subsample based on the excess 
stellar light beyond that expected from an extrapolation 
of either a smooth r 1 / 4 or an exponential functio n. Fit- 
ting a generalized exponential (or Sersic function; ISersid 
1968) to the azimuthally averaged light profiles of the gas 
indicates disklike distributions, where the threshold of a 
disk profile is taken to be n = 2 (|Cresci e t al. 2005]) and 
n = 1 is an exponential and n — 4 is a de Vaucouleurs 
profile. Table [5] lists the best-fit n values for each of the 
galaxies observed, which on average give n = 1.6±0.4. 

3.2. Kinematics and Mass 

In general, the nuclear H2 kinematics in the sam- 
ple galaxies exhibits ordered rotation with kinemetry 
fits that are consistent with co-planar circular rotation 
(i.e. higher orders of the harmonic expansion are within 
the velocity and a measurement errors). Furthermore, 
kinemetry indicates that both the PA and inclination 
angle remain constant throughout the nuclear region (to 
within the measurement error of ~10°), and residuals of 
fits with the PA and inclination angle held constant are 
consistent with the estimated kinematic errors (< 15 km 
s _1 ). This suggests that, in the galaxies observed, the 
nuclear gas has no significant radial motions (e.g. in/out 
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Fig. 1. — 2-D maps of the H2 flux distribution, velocity, and velocity dispersion, from left to right, for each Seyfert galaxy observed with 
SINFONI. In all maps of galaxies observed with SINFONI, North is alighted with the y-axis and East is to the left. The circled cross marks 
the location of the nonstellar continuum. 
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Fig. 2. — 2-D maps of the H2 flux distribution, velocity, and velocity dispersion, from left to right, for each Seyfert galaxy observed with 
OSIRIS. The y-axes in maps of NGC 3227 are aligned along a PA of -45°, and those of NGC 7469 are alighted along a PA of -48°. In all 
other maps North is aligned with the y-axis and East is to the left. The circled cross marks the location of the nonstellar continuum. 



flows) and that there is no detectable warp down to the 
smallest scales measured, typically ~20 pc. Although de- 
viations of up to ~30 km s _1 from co-planar disk rotation 
are seen in NGC 1097, NGC 3227, and NGC 4151 (see 
the Appendix for details), the kinematics in each case is 
still dominated by disk rotation. We therefore conclude, 
based on the kinematics and spatial distribution of H2, 
that the properties of the nuclear molecular gas in this 
AGN sample are consistent with a disklike distribution 



with a radial scale on the order of a few tens of parsecs. 

The best-fit PA and inclination angles given by the 
kinemetry fits to the H2 gas kinematics are listed Table 
[2] (see the Appendix for a discussion of the individual 
galaxies). In all galaxies where a comparison is possible, 
the best-fit parameters to the gas kinematics are consis- 
tent with those determined from the stellar kinematics 
measured using the CO bandhcads detected in the same 
i*T-band data (except in NGC 1068, as discussed). Values 
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Fig. 3. — Azimuthal averages based on 2-D measurements made with SINFONI. As suggested by the labels, the plots show the H2 flux 
distribution (diamonds), the Sersic fit to this distribution (dash-dotted curve), the PSF (dashed curve), the inclination corrected rotational 
velocity (solid curve with triangle data points), and the velocity dispersion (solid curve with "x" data points). The galaxy name is given 
in the upper right of each plot. 



based on the stellar kinematics for the SINFONI subsam- 
ple are those reported in D07, and those for the OSIRIS 
subsample are also determined via kinemetry and are 
reported for the first time here (the nuclear stellar prop- 
erties of this subsample will be discussed in more detail 
in a future publication). As stated, all azimuthally aver- 
aged properties are computed with the best-fit PA and 
inclination angle determined from the stellar kinematics, 
when possible. However, there is no significant difference 
in the averages computed using the best-fit values deter- 
mined from the H2 kinematics. On larger scales, others 
have also found that the orientation of the gas disk agrees 
with that of the stars. For example. iDumas et al.l (|2007l ) 
find agreement between the gas (in this case [O III] and 
H/3) and stellar disk orientation in the central kpc of a 
sample of local active and inactive galaxies. Moreover, 
the similarity of the H2 kinematics and distribution with 
that of the nuclear stellar component suggests that the 
gas and stars in the central tens of parsecs are spatially 
mixed. 

While the 2-D velocity fields of the nuclear H2 gas indi- 
cate rotating disks, the velocity dispersion suggests that 
these disks are relatively thick. On average, the rota- 



tional velocity to velocity dispersion ratio, V ro t/a (where 
V ro t—V bs/ sin i), is 0.9±0.3 at a radius of 30 pc (Fig. [7|). 
This rather low value of V ro t/a indicates that random 
motions play a significant role at these scales, despite the 
general rotation of the gas. This conclusion is drawn from 
the fact that the kinetic energy of the random motion is 
of order 3<r 2 (assuming the dispersion arises from macro- 
scopic motions and that we are observing one dimension 
of an isotropic distribution) , which is significant in com- 
parison to the energy on the order of Vj! ot from rotation. 
Consequently, when determining a dynamical mass, the 
significant velocity dispersion of the gas must be consid- 
ered. To achieve this we estimate the dynamical mass by 
Mdyn = (V^ oi + 3a 2 )r/G, where r is the radius from the 
nucleus and G is the gravitational constant. The mea- 
sured kinematics and estimated dynamical masses within 
the radii of 30pc and the H 2 light distribution HWHM 
are given in Table [3] for each galaxy, and the mean a 
values within these radii are given in Table |H The sig- 
nificance and origin of the relatively high nuclear velocity 
dispersion observed in these galaxies will be discussed in 
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Fig. 4. — Azimuthal averages based on 2-D measurements made with OSIRIS. As suggested by the labels, the plots show the H2 flux 
distribution (diamonds), the Sersic fit to this distribution (dash-dotted curve), the PSF (dashed curve), the inclination corrected rotational 
velocity (solid curve with triangle data points), and the velocity dispersion (solid curve with "x" data points). The galaxy name is given 
in the upper right of each plot. For NGC 4151, the azimuthal average of the measured flux distribution, including the bright knots of 
emission, is indicated by the open diamonds with dashed error bars, while the azimuthal average of just the extended emission is shown by 
the solid diamonds with solid error bars fsee lA,5l for details). 



troscopy with the higher spatial resolution kinematics 
measured with SINFONI indicates that the velocities 
measured in the central 0".5 are consistent with the gen- 
eral galactic rotation (Fig. EJ. This supports, in at least 
those galaxies observed with ISAAC, the suggestion that 
the nuclear gas is distributed in a rotating disk. In ad- 
dition, the consistency of the nuclear kinematics with 
the large-scale disk orientation reported in the literature 
for the sample galaxies indicates that these nuclear disks 
are well aligned with the host galaxy. In each of the 
three galaxies measured with ISAAC, the velocity dis- 
persion increases within the central few arcseconds, with 
the outer disk a(r > 5") = 43-48 km s _1 and the in- 
ner disk a(r < 0".5) = 63-96 km s _1 . These elevated 
a values are more than can be accounted for by beam 
smearing of the central velocity gradient, and therefore 
represent a real increase in a in comparison to the nor- 
mal star- forming disk measured at radii outside of 5". 
For comparison, Table [4] gives the mean a measured for 
each of these galaxies within the outer disk (r > 5") and 
inner r < 0".5 from the seeing-limited single-slit ISAAC 



data, as well as that from the inner r < 0".5 of the inte- 
gral field SINFONI data. 

For those galaxies in which the cold molecular gas has 
been measured through millimeter intcrferometric CO 
observations on scales similar to those measured here, 
the kinematics of the relatively hot H2 and cold gas are 
found to be qualitatively similar. This indicates that 
nongravitational forces do not significantly influence the 
hot gas traced by H 2 and that it does indeed follow the 
mass distribution in the nuclear region. The best-fit PA 
and inclination angles of the H2 gas disks agree with 
those of the CO gas disks reported in the literature (see 
Table [S] for references), with differences of <10 km s _1 
and thus within the measurement errors. In addition, 
similar to the hot molecular gas, the velocity dispersion 
of the cold gas is significant with respect to its rotational 
velocity (Table [5]). For example, in NGC 7469 the veloc- 
ity dispersion measured (after taking into account beam 
smearing of the velocity gradient) is as high as 60 km 
s _1 compared to a rotational velocity of ^100 km s _1 . 
Although the velocity dispersion of the cold molecular 
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Fig. 5. — Observed H2 1-0 S(l) velocity (triangles) and velocity dispersion ("x" data points) along the major (left column) and minor 
(right column) axes as measured with ISAAC. The PA of the slit (East of North) and the galaxy name are as indicated at the top of each 
plot. A comparison of the ISAAC data with the same region measured with SINFONI is shown in the smaller inset plot, where the higher 
resolution velocity and velocity dispersion are indicated by the solid curves. 



gas is generally lower than that of the hotter H2 gas, its 
relatively high value indicates that the cold gas, similar 
to the hot gas, has a significant component of random 
motion. 

3.3. Vertical Disk Height 

With a velocity dispersion that is greater than the 
inclination-corrected rotational velocity out to radii of 
tens of parsecs (Fig. [7]), a simple thin disk of H 2 gas 
can be ruled out. As discussed in H3.2[ this is because, 
in this case, the kinetic energy of the random motion is 
greater than that from the ordered rotation. In addition, 
the kinematics of the sample galaxies is consistent with 
co-planar circular rotation, suggesting that disk warps 
are not the cause of the high velocity dispersion. There- 
fore, although the molecular gas exhibits bulk rotation, 
on these scales it must be in a more spherical-like distri- 
bution (e.g. a geometrically thick disk). 

As discussed in § 13.21 the a measured in the extended 
stellar disks (r > 5") of the galaxies measured with 
ISAAC indicates that processes associated with typical 



disk star formation are capable of producing an H2 ve- 
locity dispersion of a ~ 45 km s _1 (comparable to the a 
of the cold molecular gas; Table [5]). An additional veloc- 
ity dispersion of ~ 30 km s _1 is then needed to produce 
the higher H2 nuclear velocity dispersions observed in the 
AGN at a radius of 30 pc (this is found by subtracting, in 
quadrature, the a observed in the outer galactic disk). A 
similar conclusion is reached when considering the mean 
a within a radius of 30 pc (Table I3J . The excess H2 ve- 
locity dispersion is greater than 25 km s _1 in five of the 
eight galaxies, with values ranging from 26-87 km s _1 
at a radius of 30 pc (9pc in Circinus). In the remaining 
three galaxies (NGC 3783, NGC 4051, and NGC 6814) 
the measured a is comparable to the value seen in the 
outer star forming disks. 

Models have shown that strong H2 1-0 S(l) emission 
can be produced via shocks of speeds 20-40 km s _1 
(e.g. iHollenbach fc McKedfl989l lBurton1ll992ft . How- 
ever, when the shock exceeds 40 km s" 1 the molecular 
hydrogen is dissociated, and thus emission is no longer 
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Fig. 6. — Azimuthal averages of the H2 flux distributions as 
indicated by the legend. The horizontal dashed line indicates the 
H WHM . For NGC 3227, the best fit Sersic function is shown (see 
§ IA.2I for details). Typically the standard deviation of the az- 
imuthally averaged flux is 5%. 
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Fig. 7. — Ratio of azimuthally averaged V ro t/v as a function 
of radius for each galaxy as indicated by the legend. The point 
where velocity is equal to <r is indicated by the dotted horizontal 
line. The typical standard deviation of the azimuthal average is 
less than 0.1. 

produced. Even if the excitation of H2 is dominated 
by shocks at the maximum speed of 40 km s _1 , and 
these shocks are randomly oriented, it is only possible 
to produce a ~ 40 km s _1 (including instrumental dis- 
persion). It is, however, possible to produce a higher H 2 
velocity dispersion with bow shocks, where the angle be- 
tween the shock front and the propagation direction of 
the bow is highly oblique and the resulting shock does 
not dissociate the H2 (examples in clude the Orion bu llets 
and Herbig-Haro ob ject HH99B; iTedds et alj 119991 and 
iGiannini et alj [20081 . respectively). In this scenario, the 
H2 1-0 S(l) velocity dispersion measured when integrat- 
ing over the molecular cloud corresponds to the cloud's 
bulk motion. Therefore, the high gas velocity dispersion 
measured is likely to be the result of bulk motion of the 
molecular clouds. 

Although it is not straightforward to derive the disk 
scale height from the nuclear velocity dispersions mea- 



NOTE. — PA is given as degrees East of North. For NGC 3227 and NGC 
7469 the first and second lines arc based on SINFONI and OSIRIS data, re- 
spectively. As discussed in §[3] NGC 1068 is excluded from the AGN sample. 
a For NGC 3783 and Circinus no stellar kinematics could be measured from 
the K-band data presented and instead listed here arc fits to the large scale 
K-band I Mulchacy et al. 1997) and optical 1 Freemen et al. 1977) photometry, 
respectively. In addition, for NGC 7469 the best-fit parameters based on CO 
2-1 molecular gas (assuming the gas ring at r = 2".3 is circular; Davics et al. 
12004 bl) are given on the first line instead of the SINFONI measured best-fit 
stellar parameters (which agree with those derived from the OSIRIS data 
shown on the second line). 

b See ij |A.2| for details of the Sersic fit to NGC 3227. 

° Measurement of direct HWHM radius from the observed azimuthal average, 
except for NGC 3227, where the HWHM is that of the Sersic function fit as 
described in 5 |A.2| 

sured, we attempt to get a rough estimate via two meth- 
ods. The first is to assume that the gas and stars are self 
gravitating, in which case the scale height of the gas disk 
can be estimated with the relationship z = a 2 /2irGY,, 
where G is the gravitational constant and £ is the surface 
density of the disk assuming a gas mass determined from 
the Md yn estimates derived in § 13.21 and the gas mass 
fraction of 10% estimated in the following section, i.e. 
§ 13.4.11 Strictly speaking this relationship should only 
be applied to a thin isothermal disk and consequently 
the relationship breaks down at the smallest radii where 
y rot jo is significantly less than unity (r < 20 pc). Tak- 
ing the measured a at a radius of 30 pc, where V r ot/& 
~ 1, the implied scale height is on average z — 40±6 
pc for the galaxies observed (excluding Circinus, which 
was measured out to a radius of only 9 pc where the 
disk height is 13 pc). At the HWHM radius of the light 
distribution, which is on average 23±11 pc, the average 
height is z — 33±16 pc. Estimates of the disk heights 
for the individual galaxies are given in Table 03 

The second method used to determine the height of the 
gas disk assumes vertical hydrostatic equilibrium, which 
implies V rot /cr ~ r/z a . For each galaxy, the estimated z Q 
at a radius of 30 pc and at the HWHM radii are consis- 
tent, to within 10 pc, with that derived above under the 
assumption of a thin isothermal disk. The average disk 
height given by this method is z a = 38±13 pc at a radius 
of 30 pc, again consistent with that derived above. We 
therefore estimate that the heights of the molecular gas 
disks in the observed AGNs are on the order of tens of 
parsecs at the radii considered, with typical values of z Q 
~ 30-40 pc at a radius of 30 pc. 
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TABLE 3 

Quantities Derived from the H2 1-0 S(l) Kinematics 



At r=30pc 



At HWHM Radius 



n 1 


V 

km s _1 


cr 

km s^ 1 


V/CT 


Zo 

pc 


m dyn 
10 7 Mq 


10 23 cm -2 


I' 

pc 


V 

km s — 1 


cr 

km s _1 


V/CT 


Z 

pc 


" l dyn 
10 7 Mq 


nff (Jg—U.l) 

10 23 cm" 2 


NGC 1097 


73 


52 


1.42 


30 


9.5 


4.2 


28.5 


68 


52 


1.32 


30 


8.4 


4.1 


NGC 3227 


59 


98 


0.60 


45 


22.6 


10.0 


29.3 


57 


98 


0.58 


41 


22.0 


10.2 




71 


82 


0.87 


40 


17.7 


7.8 


34.9 


86 


87 


0.99 


41 


24.3 


7.9 


NGC 3783 


26 


33 


0.78 


40 


2.8 


1.3 


20.2 


6 


34 


0.18 


33 


1.6 


1.6 


NGC 4051 


37 


44 


0.83 


41 


5.0 


2.2 


7.6 


13 


49 


0.26 


13 


1.3 


8.9 


NGC 4151 


82 


67 


1.23 


33 


14.0 


6.2 


28.1 


76 


59 


1.28 


30 


10.8 


5.1 


NGC 6814 


23 


13 


0.53 


46 


4.2 


1.9 


34.8 


26 


47 


0.55 


57 


6.3 


1.8 


NGC 7469 


34 


64 


0.54 


46 


9.3 


4.2 


77.4 


68 


71 


0.96 


99 


35.4 


2.4 




43 


62 


0.69 


43 


10.1 


4.2 


30.8 


45 


63 


0.71 


43 


10.7 


4.2 


Circinus a 


36 


51 


0.70 


13 


1.9 


9.4 


7.4 


38 


54 


0.71 


11 


1.8 


12.9 


NGC 1068 


21 


102 


0.21 


49 


22.5 


10.1 


13.2 


8 


118 


0.07 


22 


12.8 


29.2 



Note. — The typical uncertainty on the velocity and velocity dispersion measurements is 5-15 km s - . For NGC 3227 and NGC 7469 the first 
and second lines are based on SINFONI and OSIRIS data, respectively As discussed in § [3] NGC 1068 is excluded from the AGN sample. The 
disk scale height, z , is derived as described in § 13.31 using the relationship z — cr 2 /27rGS. Mdyn and Nh are derived assuming the typical 
10% gas mass fraction - sec §3.2-3.4 for details and discussion. 

a Quantities measured at 9 pc (the edge of the measured FOV) rather than 30 pc. 



TABLE 4 

Mean H2 1-0 S(l) Velocity Dispersion 



SINFONI & OSIRIS 



ISAAC 



Galaxy 


<<T>30pc 


<<J>HWHM 


<(T> If 

.0 


r 


<0" >D07b 


Cmax 


<cr> 2o" 


<cr> // 


<a> 


".5 




km s _1 


km s 


km s _1 


pc 


km s — 1 


km s — 1 


km s 


km s 


km s 


-1 


NGC 1097 


61± 9 


62±10 


60± 9 


22 


64± 9 


74 


58±10 


47±11 


71± 


1 


NGC 3227 


95±15 


95±15 


94±14 


32 


95±15 


108 


74±15 


43±23 


96± 


5 




78±14 


79±14 


77±14 




79±14 


93 










NGC 3783 


33± 1 


33± 1 


34± 2 


60 


33± 1 


33 










NGC 4051 


55± 8 


56±12 


56± 8 
















NGC 4151 


56±10 


54±10 


57±11 
















NGC 6814 


48± 6 


48± 5 


45± 6 
















NGC 7469 


62± 2 


67± 5 


71± 7 


128 


70± 7 


82 












62± 6 


62± 6 


79±12 




78±13 


92 










Circinus a 


56± 3 


57± 2 


56± 3 


8 


56± 2 


61 


54± 9 


48±10 


63± 


4 


NGC 1068 


88±42 


58±40 


88±42 


30 


88±42 


134 











Note. — A uniformly weighted mean within a radius given by the subscript, except where it is indicated to be outside of the 
stated radius. For NGC 3227 and NGC 7469 the first and second lines are based on SINFONI and OSIRIS data, respectively. 
As discussed in §[3] NGC 1068 is excluded from the AGN sample. 
a Quantities measured at 9 pc (the edge of the measured FOV) rather than 30 pc. 



TABLE 5 

Summary of CO Measurements 



Galaxy 


Radius 


Mgas 


Mdyn 


h 


cr 


v/cr 


References 




tt 


PC 


10 7 Mq 


10 7 Mq 


% 


km s _1 






NGC 1097 


1.8 


158 


2.6 


207 


1.3±0.2 


50 


5 


1 


NGC 3227 


0.3 


24 


0.8±0.3 


6.4 


13±4 


60 


0.8 


2 


NGC 7469 


2.5 


800 


380 


650 


58 


60 


0.6 


3 


NGC 1068 


0.5 


35 


2.3-4.0 


13.6 


17-29 


30 


1.2 


1 



Note. — R eferences for CO m eas urements and derived M oaB and Mdyn are 
as foll ows: fl) iHsieh et al.l|2008l. (2) iSchinnerer et al.ll2000al. (3) iDavies et all 
20043, f4) iSchinnerer et al.ll2000bl . iMuller Sanchez et al.ll2009l . D07b. 
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TABLE 6 

Gas Mass Estimates as Derived in ^3.4.11 




/ 9 =10% Lh2 Conversion a 



Galaxy 




M sas 


M sas 


fg 




10 7 M 


10 7 Mq 


10 7 M 


% 


NGC 1097 


9.5 


1.0 


0.5-4.5 


5-47 


NGC 3227 


22.6 


2.3 


2.4-22.2 


11-98 




17.7 


1.8 


2.2-20.2 


>12 


NGC 3783 


2.8 


0.3 


1.1-10.6 


>39 


NGC 4051 


5.0 


0.5 


0.7-6.6 


>14 


NGC 4151 


14.0 


1.4 


1.2-10.8 


9-77 


NGC 6814 


4.2 


0.4 


0.6-5.5 


>14 


NGC 7469 


9.3 


0.9 


4.9-45.5 


>53 




10.1 


1.0 


4.7-43.3 


>47 


Circinus 


1.9 


0.2 


0.2-1.5 


11-79 


NGC 1068 


22.5 


2.3 


4.2-39.3 


>19 



Note. — Dynamical masses are based on the 
H2 kinematics as discussed in i)3.2l All masses are 
derived for a radius of 30 pc, with the exception 
of Circinus, for which the radius is 9 pc. For 
NGC 3227 and NGC 7469 the first and second 
lines are based on SINFONI and OSIRIS data, 
respectively. As discussed in § [3] NGC 1068 is 
excluded from the AGN sample. 
a For cases where the typical conversion from H2 
luminosity to M 9as results in a mass higher than 
Mdyn, only the f g lower limit derived from the 
more extreme conversion factor is given (see the 
discussion in £|3.4.1[1 . 

3.4. Column Density 
3.4.1. Column Density Derived from Dynamical Mass 

An order of magnitude estimate of the gas mass, and 
thus the column density, within the nuclear regions of the 
galaxies can be derived by assuming that the nuclear gas 
masses are a fixed fraction of the dynamical masses. In 
addition, we make the simplifying assumption that the 
H2 gas mass at a given radius is distributed uniformly. 
The derived column densities are thus the average val- 
ues within a given a radius, and provide a measure of 
the maximum extinction. As will be shown in the fol- 
lowing section, a significant fraction of the gas within 
the thick disks observed is in a clumpy distribution and 
the actual column densities are consequently lower than 
derived here. 

This method of estimating the column densities re- 
quires that the gas masses are distributed as implied 
by H2 1-0 S(l) and are not in a denser component 
such as the H2O mas er disks seen on parsec scales (e.g. 
IGreenhill et all[200l . We argue that this is indeed the 
case based on evidence that the nuclear stellar compo- 
nent is spatially mixed with the thick gas disk observed 
(§ 13. 2p . The thick stellar disk indicates that the cold 
component of gas from which the stars formed must also 
be in a similarly thick disk, which is supported by the 
high velocity dispersion measured in CO observa tions 
(e.g. ISchinnerer et all l2000al iDavies et all l2004ri Ta- 
ble Therefore, while evidence exists of dense thin 
disks on scales smaller than observed here, the distribu- 
tion and kinematics of the stars and molecular gas (both 
hot and cold) suggests that, on scales of tens of parsecs, 
the mass is distributed in a thick disk consistent with the 
disk traced by H 2 1-0 S(l). 



m a l 

20 40 60 80 100 120 140 
radius (pc) 

Fig. 8. — Azimuthal average of the average gas column density 
for each galaxy as indicated by the legend. The average column 
density is calculated assuming that the gas is uniformly distributed 
throug hout t he region considered and that t he gas fraction is 10% 
(see § 13.4.11 for details). As discussed in § 13.4.11 the molecular 
gas is actually in a clumpy distribution, and the densities plotted 
here represent only the averaged column densities within the given 
radius. The error bar shows the typical 35% error of the gas column 
density estimates. 

Although the gas mass fraction (f g ) is not accurately 
known in the nuclear regions of AGN, we attempt to 
constrain this value via several methods discussed be- 
low. I. Typical Galaxies — Normal spiral galaxies, sim- 
ilar to the galaxies hosting the Seyfert AGN observed, 
have typical gas mass fractions of ~10%, with values 
ranging from 4% for S ab to 25% for Scd galaxies (e.g. 
I Young fc Scovilld li~991). Given the evidence reported by 
D07 for recent nuclear star formation in the SINFONI 
galaxy subsample, the gas mass fraction measured in 
starburst and ULIRG galaxies may also be applicable to 
the nuclear regions of these AGNs. Local starburst and 
ULIR G galaxies are typically found to have f„ = 10-20% 
(e.g. I Young fc Scovilld fl98l IDownes fc Solomon! fl998l 
Genzel et all I2001L iGreve et al.l 120051 iTacconi et al.l 
2006). Therefore, a gas mass fraction similar to that in 
spiral, starburst, or ULIRG galaxies suggests f g ~ 10% 
with a range from 4% up to 20-25%. 
2. Cold Molecular Gas — Measurements of the cold 
molecular gas within the nuclear region exist in the lit- 
erature for four of the galaxies in this sample: NGC 
1097, NGC 3227, NGC 1068, and NGC 7469. The 
standard galactic CO-to-H2 conversion factor is not as- 
sumed since evidence indicates that this conversion fac- 
tor is lower in intense star for ming environments such as 
the nuclei observed here (e.g . IDownes fc Solomonlll998l 
iPapadopoulos fc~A llcn 2000). We therefore use the con- 
version f actor of (1.38±0.46)x l0 20 cm" 2 (K km s" 1 )" 1 
found by IDavies et alj (|2004bD for the central 800 pc in 
NGC 7469, which is 0.4-0.8 times the Galactic conversion 
measured bv lStrong et al.l ([1988). The total gas mass is 
then determined from the H2 gas mass implied by the 
CO luminosity by assuming a mass correction for helium 
of 40% and that the masses of H I and H II are negli- 
gible at these radii. As can be seen in Table O the gas 
fractions derived from the CO 2-1 integrated intensity 
vary significantly from source to source, with values of 
f g = 1-29%. To date, there are no higher spatial reso- 
lution CO 2-1 observations of NGC 7469 than those of 
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iDavies et ail (|2004bf) . but over this relatively large region 
of 1.6 kpc, the independently determined H2 mass (used 
to estimate the CO-to-H2 conversion factor above) and 
the CO dynamical mass imply f g — 58%. 
3. Sample Wide Average for Active Galaxies — There 
now exist in the literature several studies of the cold 
molecular gas in the central hundreds of parsecs of galax- 
ies of various activity levels, including Seyfert galax- 
ies. Taking a sample of five active galaxies that have 
been observed in CO 2-1 as part of the NUGA Sur- 
vey (NGC 4321, NGC 4579, NGC 4826, NGC 6574, 
and NGC 6951) we use the reported gas and dynamical 
masses to derive a gas mas s fractio n dLindt-Krieer et ah 
20081 iGarria-Burillo et all l2005l IGarria-Burillo et al. 
20031 iKrips et al. 2007). In several cases no dynamical 
mass is given in the literature at the small radii consid- 
ered here, in which case it is estimated from the reported 
CO 2-1 kinematics. Based on this sample we find a typ- 
ical gas mass fraction of ~10%, with f g = 1-15% within 
the central ~200 pc of four galaxies and f g — 10% out to 
a radius of 300 pc in the fifth galaxy. 
4. H 2 1-0 S(l) Luminosity — The gas mass can also be de- 
termined (to within a la uncertainty of a factor of about 
2) from the lumi nosity of H2 1-0 S(l) using t he relation- 
ship reported bv lMuller Sanchez et al.1 (|2006l ) , which has 
a con version factor of 4000 M0/L0 (also see lDale et al.1 
2005). A lower limit on the gas mass is estimated from 
the most extreme case in their sample, NGC 6240, which 
is known to be over luminous in H2 1-0 S(l) (e.g. 
ISugai et"al"1ll997llTecza et al.ll2000h . and has an order of 
magnitude lower conversion factor of 430 M /L Q . Based 
on this extreme conversion factor, the median gas mass 
fraction for the sample is f g ~ 13%, with values rang- 
ing from 5-53% (individual gas mass estimates are given 
in Table O- In each of the individual galaxies, the gas 
mass fraction derived from the H2 luminosity is consis- 
tent with that determined from the CO observations dis- 
cussed above. The typical conversion factor gives sig- 
nificantly higher gas masses and the lowest value in the 
sample is f g — 47%. In some cases the typical conversion 
factor implies a gas mass fraction greater than 100%, 
suggesting that these galaxies may also be over lumi- 
no us in H2 compared to the sa mple of galaxies studied 
bv iMuller Sanchez et all (|2006h . and that a conversion 
similar to that found for NGC 6240 is more appropriate. 
We therefore deduce from this method that f g is greater 
than a few percent, with no constraint on the upper limit. 

To summarize, typical galaxies have f g ~ 10% (ranging 
from 4 to 25%), CO measurements of active galaxies also 
suggest f g ~ 10% (with a range from 1 to 15%), CO mea- 
surements of three Seyferts in this sample on scales com- 
parable to those measured here imply f g — 1-29%, and a 
conversion of the nuclear H2 luminosity observed in the 
sample galaxies gives a lower limit on f g of a few percent. 
We therefore conclude that for the nuclear regions of the 
Seyfert galaxies observed that typically f g ~ 10% with a 
range of f g ~ 1-20%. If the gas is uniformly distributed, 
then f g ~ 10% gives an average column density within 
a radius of 30 pc of N H ~ (4.9±3.3) x 10 23 cm" 2 for 
the galaxy sample, with values ranging from 1.3 to 10.0 
x 10 23 cm -2 (Fig. [5]). The column densities within a 
radius of 30pc and the HWHM light radius are listed for 
each galaxy in Table [3] Taking f g ~ 1% still results in 



column density estimates greater than 10 22 cm 2 in all 
galaxies, which is sufficient to obscure an AGN, and f g ~ 
20% gives an average Nh ~ 10 24 cm~ 2 . However, as will 
be shown in the following section, the gas is unlikely to 
be uniformly distributed within the nuclear region, and 
the actual clumpy distribution results in significantly less 
obscuration of the AGN (as well as the nuclear star for- 
mation) . 

3.4.2. Column Density Derived from Extinction 

The column densities within the nuclear regions of the 
Seyfert galaxies can also be determined from the extinc- 
tion of the nuclear stars, which can be constrained by 
fitting a reddened stellar template to the nuclear stellar 
continuum. We use K-type stellar templates, but the 
template used does not significantly alter the derived 
column densities because in the if-band there is little 
difference in the intrinsic continuum slope for any rea- 
sonable stellar population. The best-fit extinctions for 
the individual galaxies are given in Table [7] assuming a 
model in which the gas and dust are mixed with the stars 
{Iobs—Io{l-e~ Tx )/T\, where I and I b s are the intrinsic 
and observed luminosities, and t\ is the optical depth 
at a given wavelength) . Extinction estimates range from 
Ay = 1 mag to 24 mag, with an average extinction of 
12 mag for the sample. This amount of extinction, as- 
suming Nh— 1.9 x 10 21 cm~ 2 mag^ 1 Ay, implies an 
average column density of Nh = (2.4±1.6) x 10 22 cm -2 , 
with values ranging from 0.2 to 4.5 x 10 22 cm~ 2 (see 
Table [7]). Also given in Table [7| is the best-fit extinctions 
and implied column densities assuming screen extinction 
(I b s =I c~ Tx ). This model provides a lower limit since, 
as discussed in § 13 . 21 the similarity of the kinematics and 
distribution of the stars and gas in these galaxies sug- 
gests they are mixed. As a result a screen model is not 
expected to be applicable to the nuclear regions of these 
galaxies. 

The nuclear stellar continuum in the Seyfert galaxies 
is extracted from an annulus with inner and outer radii 
of 25 and 35 pc, and the region measured is therefore 
matched to the r ~ 30 pc spatial scales of the nuclear 
gas disks. Although the high spatial resolutions achieved 
with the AO typically gave PSFs with an HWHM of 
< 0".15 (less than 16 pc at the average distance of the 
galaxies observed), in some galaxies the contamination 
within this annulus from the nonstellar AGN continuum 
in the PSF wings is nonnegligible. To account for this 
contamination, the nuclear spectrum extracted from a 
region dominated by nonstellar continuum (within the 
FWHM of the PSF core; Table [lj is first subtracted from 
the annular spectrum. The nuclear spectrum is scaled to 
match the annulus spectrum broad Br7 emission, which 
is assumed to originate in the unresolved BLR. The com- 
bination of the formal uncertainties given in Table [7| and 
the uncertainty introduced by error in the correction for 
the AGN continuum contamination, as well as template 
mismatch, results in an overall accuracy of the extinction 
estimates of up to ~50% with typical values of ^30%. 
Correction of the continuum for nonstellar contamina- 
tion is highly uncertain in the two Seyfert 2 galaxies be- 
cause there is no broad Br7 emission to use for scaling 
of the nuclear spectrum. A correction is attempted by 
scaling the nuclear spectrum by the flux level in the PSF 
wings at the radii measured, and these estimates are re- 
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TABLE 7 
Best-Fit Extinction 



Screen Mixed 



Galaxy 


r 






A . 






pc 


mags 


10 21 cm -2 


mags 


10 21 cm -2 


NGC 1097 


25-35 


3.8 


7.2 


8.8 


16.7 






(0.5) 


(1.0) 


(1.5) 


(2.8) 


NGC 3227 


25-35 


4.1 


7.7 


9.7 


18.5 






(0.3) 


(0.7) 


(1.0) 


(1.9) 




25-35 


4.3 


8.1 


10.2 


19.3 






(0.4) 


(0.7) 


(1.1) 


(2.0) 


NGC 3783 


60-80 


8.5 


16.1 


27.7 


52.7 






(0.1) 


(0.2) 


(0.6) 


(1.1) 


NGC 4051 


25-35 


0.0 


0.0 


1.0 


1.8 






(0.4) 


(0.7) 


(1.1) 


(2.1) 


NGC 4151 


25-35 


5.6 


10.7 


14.5 


27.6 






(0.6) 


(1.2) 


(2.2) 


(4.2) 


NGC 6814 


25-35 


1.4 


2.7 


3.0 


5.7 






(0.4) 


(0.8) 


(0.9) 


(1.7) 


NGC 7469 


60-80 


9.5 


18.0 


37.1 


70.5 






(0.1) 


(0.2) 


(1.4) 


(2.7) 




60-80 


6.5 


12.4 


17.7 


33.6 






(0.2) 


(0.4) 


(0.8) 


(1.5) 


Circinus 


7-9 


7.5 


14.2 


22.5 


42.7 






(0.3) 


(0.5) 


(1.6) 


(3.1) 


NGC 1068 


25-35 


12.8 


24.3 


121+ 


230 1 " 






(0.4) 


(0.7) 


( ••■ ) 


( ■■• ) 



Note. — Values in parenthesis are the error of the measurement 
immediately above. For NGC 3227 and NGC 7469 the first and 
second set of lines are based on SINFONI and OSIRIS data, re- 
spectively. For the Seyfert 2 galaxies Circinus and NGC 1068 the 
correction for potential non-stellar contamination is uncertain (see 
text for details). The values marked with f for NGC 1068 are be- 
yond the saturation limit of the mixed extinction model, which is at 
an optical depth of a few. As discussed in §[31 NGC 1068 is excluded 
from the AGN sample. 

ported in Table [7] along with the uncorrected estimates. 
The corrected value for Circinus is used when determin- 
ing the sample-wide properties, and, as discussed, NGC 
1068 is excluded from the sample. Therefore, the general 
conclusions of the stellar extinction in the sample are not 
affected by the additional uncertainty of the corrected es- 
timates in the Seyfert 2 galaxies. In addition, due to the 
greater distance of NGC 3783 and NGC 7469, an esti- 
mate of the extinction from r = 25-35 pc was not pos- 
sible because the PSF HWHM already extends in both 
cases to >20 pc resulting in significant nonstellar con- 
tamination of the continuum at r ~ 30 pc. Instead, the 
extinction is measured for these two galaxies in a larger 
annulus of r — 60-80 pc. Therefore, for these two galax- 
ies, the estimates are likely to represent a lower limit on 
the extinction at r ~ 30 pc. Also, for Circinus the ex- 
tinction within r = 7-9 pc is given in Table [7l which is 
at the edge of the FOV. 

The relatively low extinction of the stellar continuum 
in these galaxies implies column densities that are over 
an order of magnitude lower on average than those de- 
rived within a similar spatial region in § 13.4.11 under the 
assumption that the gas comprising a fixed fraction of 
the dynamical mass and is uniformly distributed (me- 
dian [extinction] = 0.04A^[H 2 / s =10%]). One way 
to reconcile this discrepancy is to assume a lower gas 
mass fraction. However in all but one case (NGC 3783) 
a value of f g < 1% is required to bring the column 
densities into agreement. This value is lower than ob- 



served in any known Seyfert or typical spiral or star- 
burst galaxy (see § 13.4. 1|) . A more likely scenario is 
that the gas within this nuclear region is in a clumpy 
distribution such that the nuclear stellar and central 
AGN emission is not heavily obscured despite the large 
gas masses present. As will be discussed in more de- 
tail in the following section (and was mentioned in § [1} , 
this clumpy distribution is also supported by theoretical 
mod els (e.g. lElitzurll2006bf) . h vdrodvnamical simulations 
(e.g. ISchartmann et aLll2009f). and observat ions of ther- 
mal dust emissi on (e.g. iTristram et al.ll2007l ) and silicate 
absorption (e.g. iDeo et al.ll2007l ). 

3.5. Relatively Transparent Thick Disks of Molecular 

Gas 

In summary, based on the kinematics of H 2 1-0 S(l), 
the molecular gas within the central ~100 pc of the 
Seyferts is in a rotating, thick (z /r ~ 1) disklike dis- 
tribution that is spatially mixed with the nuclear stellar 
population. This gas/stellar disk has spatial scales of less 
than 35 pc in all cases, with a typical HWHM of ^30 pc. 
Within the central ~60 pc the dynamical mass is ^10 8 
M Q , which suggests M gas ~ 10 7 M Q based on the typ- 
ical gas mass fraction of 10% (M gos = (0.1-2.0) x 10 7 
Mq assuming f g = 1-20%). If a uniform gas distribution 
is assumed then this typical gas mass fraction implies an 
average column density at a radius of 30 pc of ~ 5 x 10 23 
cm~ 2 . However, the average column density implied by 
the extinction of the nuclear stellar continua is over an 
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order of magnitude lower than this estimate (Njj ~ 2 
x 10 22 cm" 2 ), suggesting that the nuclear gas in these 
galaxies is actually in a clumpy distribution. The disklike 
distribution of gas on — 10 pc scales in Seyfert galaxies is 
therefore relatively transparent and does not heavily ob- 
scure the nuclear stellar population or the central AGN. 

In addition to theoretical models and hydrodynam- 
ical simulations indicating that the gas on ~10 pc 
scales is distributed in discrete clouds rat her than evenly 
throughout the nuclea r region (e.g. lElitzuri [2006b, 
iSchartmann et aT]|2009f ). additional support for this con- 
clusion comes from spectroscopic observations of 10 /im 
silicate in Seyfert galaxies. It is understood that signifi- 
cant silicate absorption requires the radiation source be 
deeply embedded in dust that is both optic ally and geo- 
metrically thick (e.g. iLevenson et al.ll2007f) . Strong sili- 
cate absorption therefore suggests that a fairly uniform 
thick shell of gas and dust surrounds a compact star- 
burst or AGN. On the other hand a medium that consists 
of clouds illuminated from the outside will not produce 
silicate absorption since there is no strong temperature 
gradient across the clo uds leading to a rela tively con- 
stant dust temperature. lHao et aiT(|2007n and lDeo et al.1 
(|2007f ) found that only a small fraction of Seyferts have 
deep silicate absorption (seven out of 61 Seyfert galaxies, 
six of which are type 2s) . The lack of silica t e absorption 
in most Seyfert galaxies dShi et aTl I2006L iSpoon et al.1 
120071 lUddin et al.ll2007l iThompson et al.ll2007l ) is inter- 
pret ed to mean that the gas and dust is clumpy. More- 
over, [DejZeFlil] (|2007h find that the few Seyfert with deep 
silicate absorption arc cither highly inclined or merging 
systems, suggesting that the spectral feature arises in a 
c ool dusty reg i on on scales of tens of parsecs or greater. 

IMarr et"aTT(|1993h find that in our own Galaxy the 
molecular gas in the central tens of parsecs is confined 
to virial clumps with size s of order 0.1 pc and den sities 
of 10 6 cm" 3 . Similarly, IChristopher et al.1 (|2005f ) find 
numerous dense cores with typical diameters of 0.25 pc, 
densities of 3-4 x 10 7 cm" 3 , and masses of ~ 2 x 10 4 
M Q . With such high densities, these clouds are stable 
against tidal disruption and have long lifetimes (~10 
yrs). In addition to these virial estimates, similar cloud 
densities are also der ived from multitransition modeling 
ijJackson et al.l [1993T ) and are also found for clouds in 
the central 30 pc of Seyferts usin g analytical modeling 
(|Vollmer. Beckert. fc Daviesl I2008D . If we assume that 
all of the molecular gas in the observed Seyfert galaxies 
is composed of clouds with properties similar to those 
seen in the Galactic Center (0.25 pc diameter, ~ 2 x 
10 4 M ), then there must be about 500 clouds within a 
radius of 30 pc for a gas mass of - 10 7 M Q . This results 
in only a ~1% chance that a given line of sight will be 
blocked by a cloud, and considering the full range of f g 
values of 1-20%, the chance of obscuration by a cloud 
is 0.1-2%. These results are qualitatively in agreement 
with the small fraction of Seyferts that show deep silicate 
absorption. However, this is an extreme case in which all 
of the gas is assumed to be within dense clumps. Since 
the reddened stellar continuum in these galaxies implies 
significant extinction, it is likely that in addition to dense 
clouds of gas there is a nonnegligible fraction of the nu- 
clear gas in a diffuse component. 



80 


• Mean a 




. i ( 


NGC 322? ' 

1 A 












A Max a 








60 






NGC 7469 . 








i 1 


> A 




40 


ircmu> 






- 


20 


TNGC 3783 , "V 

y i 


NGC 1097 

• A 



















40 


60 

a (km s~ 


80 
') 


100 



Fig. 9. — H2 velocity dispersion vs. the density of the time- 
averaged SFR (averaged over the time elapsed since the beginning 
of the star-forming episode), as determined by D07 and adjusted 
so all galaxies are averaged over 200 Myr. For each galaxy, a circle 
shows the mean dispersion within the radius measured by D07 and 
a triangle to the right of each data point indicates the maximum 
dispersion within this radius. The horizontal error bars show the 
standard deviation of the velocity dispersion and does not take into 
account the typically 5-15 km s — 1 measurement errors. Fits to 
the mean and max velocity dispersions with respect to the time- 
averaged SFRs are shown by the dashed and dash-dotted lines, 
respectively. NGC 1068, which is excluded from the sample (see 
§E), has a time-averaged SFR of 163 M Q kpc~ 2 over 200 Myr, 
placing it outside of the region plotted. 

4. CORRELATION OF THE MOLECULAR GAS AND STAR 
FORMATION PROPERTIES 

We find a correlation between the H2 1-0 S(l) velocity 
dispersion and the SFR density in the sense that galaxies 
with greater star formation have higher velocity disper- 
sion, implying a thicker gas disk. This relationship can 
be seen in Fig. [9j where the time-averaged SFRs deter- 
mined by D07 are all normalized to 200 Myr. Plotted for 
each galaxy are both the mean and maximum a of H2 
within the region considered by D07 (0".25-0".4, or 8-128 
pc, depending on the galaxy; see Table H]). Linear fits to 
both of these quantities with respect to the SFR den- 
sity indicate that there is a significant correlation, with 
linear correlation coefficients of —0.8. The significance 
of this correlation and the suggested link between star 
formation and height of the gas disk will be investigated 
further in the following section, i.e. § [5l 

A correlation is also found between the surface density 
of the molecular gas and the SFR per unit area that 
is consistent with known scaling relations. As shown 
in Fig. [TTJ the observed AGN for which SFRs have 
been determined generally agree with both the Schmidt- 
Kennicutt relation, which assumes a co nstant CO-to-Hg 
conversion typical of normal galaxies (|Kennicuttl I1998L 
Kennicutt et al.ll2007f) . and with the modified relation of 
Bouche et al.ll2007l which accounts for the decreased CO- 
to-H2 conversion found in denser and more active sys- 
tems. Based on the time-averaged SFRs (D07) shown in 
Fig. EH the AGN are in better agreement with the modi- 
fied relation. Although the current SFRs, which are typ- 
ically considered in this relation, are not well constrained 
in the observed AGNs, they are likely to be at least an 
order of magnitude lower than the time-averaged SFRs 
(D07). Applying this correction, the AGNs are then 
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Fig. 10. — Gas surface density vs. the SFR per unit area (aver- 
aged over the time since the beginning of the star-forming episode, 
see D07 for details) for the SINFON I subsample galaxies co mpared 
to the Schmidt-Kcnnicutt relation (Kcnni cutt et all 1200 7; dashed 
line with the formal error of the fi t indicated by dotted lines) and 
the adjusted relation reported by Bouchc ct al. (2007) (solid line 
with the fit error indicated by the dash-dotted lines) . The gas sur- 
face density at the radii measured by D07 are shown by the circles, 
and to the right of each data point the mean surface density within 
this radius is indicated. The data shown assume f g = 10%. If f g = 
1% is assumed, then all galaxies have gas surface densities that are 
lower at a given SFR per unit area than either relation. The scat- 
ter of the galaxies observed is similar to that se en in other samples 
(see, for example, Fig. 3 of Bouchc ct al. 2007). 

found to lie closest to the orig inal Schmidt-Kennicutt re- 
lation of lKennicutt et aT]|2007l . Furthermore, a gas mass 
fraction of 10% is assumed to derive the gas surface den- 
sities shown in Fig. [TUJ If a gas mass fraction of only 1% 
is assumed, which, as discussed in § 13.4.11 is also a plau- 
sible value for the nuclear regions of these galaxies, then 
all galaxies have a lower gas surface density for a given 
SFR per unit area than either relation. The correla- 
tion of the gas surface density with the SFR density, and 
general agreement with the relations observed in other 
galaxies, suggests that the high SFRs observed in these 
AGNs should not come as a surprise given the high gas 
surface densities in the nuclear regions. 

5. ACHIEVING THE IMPLIED DISK HEIGHT 



As discussed in § 13.31 the high velocity dispersion with 
respect to the rotational velocity implies a geometrically 
thick gas disk. However, the individual molecular clouds 
of the thick disk would cool via supersonic shocks from 
collisions on a relatively short timescale (comparable to 
the orbital period of a few million years), thus energy 
must continually be injected into the region in order 
to maintain this vertical structure. Many theories have 
been proposed in the literature to address how the verti- 
cal structure of the obscuring medium, or "torus" , can be 
maintained. The feasibility of these mechanisms in main- 
taining the vertical structure on scales of tens of parsecs 
in the observed Seyfert galaxies is now considered. 

As discussed in § 13.21 kinematically the molecular hy- 
drogen in the galaxies observed is consistent with a ro- 
tating disk and there is no evidence of a significant warp 
(e.g. twisting of the kinematics major axis PA and/or 
inclination angle) in the galaxies down to scales of ~20 
pc, nor of significant radial motion such as an outflow. 



Therefore, for the galaxies in this sam ple, models of disk 
winds (e.g. lElitzur fc Shlosmanl 20061) and warped disks 
(|Navakshinl 120051 ICaproni et al. 120061 ) are ruled out as 
the dominant source of the high velocity dispersion at 
the radii considered here, and are thus not viable mech- 
anisms to explain the implied vertical height of the gas 
disk. It has, however, been suggested that there is evi- 
denc e of warped gas disks on these sca l es in other AGN 
(e.g. iKrajnovic. Sharp, fc Thattdl2007l iNeumaver et all 
120071) . and, based on H 2 Q maser observations, on smaller 
< 1 p c scales (e.g. iHerrnstein et al.| [l996. Grc enhill et al.l 
2003). In addition, a disk wind may still be present on 
these scales but if the surface brightness is too low then it 
would be undetected. In this case, the disk wind would 
not be a significant component of the H 2 emission de- 
tected and would thus not appreciably impact the H 2 
distribution nor kinematics measured. 

Based on stellar continuum intensity of the galaxies 
also measured with ISAAC, the nuclear starburst has a 
significantly higher SFR in comparison to the star for- 
mation of the outer disk (r > 5"), where a velocity dis- 
persion of only ~45 km s _1 is measured (see § 13 . 2[) . The 
contribution from star formation to the velocity disper- 
sion, and therefore to the disk height, is thus likely to 
be greater in the nuclear region, especially since the sim- 
ilarity of the distribution and kinematics of the molecular 
gas and stars suggest they are spatially mixed. The nu- 
clear star formation can contribute to the H 2 disk height 
through at least three different mechanisms: stellar ra- 
diation pressure, supernovae, and stellar winds. 

The model of Thompson et al. (2005; see also 
iBallantvne] l2008f) suggests that stellar radiation pres- 
sure can result in disk heights comparable to the ra- 
dial distance, provided the SFR is sufficiently high. In 
their optically thin model SFR scales as SFR ~ fg a s a4 ' 
and in their optically thick model SFR ~ f ga s& 2 , where 
fgas is the gas fraction and a is the gas velocity disper- 
sion. Assuming a gas fraction of 10% (see § 13 .4. 1 [) . the 
time-averaged SFRs reported in D07 are high enough to 
achieve the disk heights implied by the a at 30 pc in both 
of these models. For example, in NGC 3227 the velocity 
dispersion of 98 km s -1 requires a SFR of at least 0.25 
Mq yr _1 to be maintained via radiation pressure from 
star formation (based on the optically thin model), and 
the time-averaged star formation is estimated to be 1.2 
Mq yr , almost five times higher than required. Con- 
sidering that the SFRs at their peak were at least a factor 
of a few up to as much as an order of magnitude greater 
than the time-averaged SFRs (D07), it is feasible, based 
on these models, that the disk height is generated via 
stellar radiation pressure. If lower gas mass fractions 
are considered then the required SFRs decrease and this 
mechanism becomes even more likely to play a role in 
maintaining the disk height. However, although D07 find 
that the nuclear star formation is recent, it is no longer 
active, and depending on the evolution of the star for- 
mation with time (which is not known) the current star 
formation could be as much as an order of magnitude 
lower than the time-averaged SFR (D07). Therefore, al- 
though this mechanism may have played a critical role 
in generating the disk height during the peak of the star 
formation and for some time afterward, it is unclear if the 
current SFRs can maintain the estimated disk heights. 
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Another possible way by which star formation can con- 
tribute to the heig ht of the gas disk i s with energy input 
from supernovae. IWada fc Norman! (|2002ft have shown 
that with turbulence due to supernovae it is possible to 
achieve disk heights of tens of parsecs. However, the 
supernovae rate required to reach these heights is ~1 
yr" 1 , which is at least one, and as much as four, orders 
of magnitude higher than the supernovae rate predicted 
for the nuclear starbursts in the SINFONI subsample 
(D07). The last mechanism related to star formation 
considered for explaining the vertical structure of the 
gas disk is stellar winds. iNavakshin fe Cuadral (2007) 
model the obscuring properties of stellar winds from a 
starburst within an accretion disk and find that a disk 
height of only a few parsecs can be achieved even with a 
significant outflow rate of 1-100 M Q yr" 1 . Although in 
principle both of these mechanisms, supernovae and stel- 
lar winds, are likely to play some role in maintaining the 
H2 disk height, neither method significantly contributes 
to the disk height of tens of parsecs implied by the high 
velocity dispersion. 

An alternative energy source capable of supporting a 
geometrically thick gas disk is radiation pressure from 
the AGN emission. Pier fc Kroliki (|1992f) show qualita- 
tively that disk heights of z o =l-10ri„ can occur when 
Lbh / LEdd=0. 01-0.1, where ri n is the inner radius of the 
torus as determined by the dust sublimation tempera- 
ture. In the galaxies observed r in ~ 0.2 pc, implying that 
this mechanism can only achieve disk heights on the or- 
der of parsecs, with the greatest height of 4 pc achieved 
in NGC 7469 . The results of the more detailed study of 
iKrolikl (|2007f ) also suggest that for the observed galaxies 
the maximum disk height is on the order of parsecs, and 
that at a radius of ~30 pc there is no significant disk 
height achieved via AGN radiation pressure. Therefore, 
although AGN radiation pressure will inevitably play a 
role in the vertical structure of the disk, it is not enough 
on its own to explain the height of the disk observed at 
the radii considered. 

An alternative to the models discussed above is pro- 
duction of the disk height through energy input from 
mass accretion int o the nuclear region. This mechan ism 
has been shown bv lVollmer. Beckert. fc Daviesl (|2008h to 
be capable of producing disk heights on the order of tens 
of parsecs with predicted inflow rates of 0.07-4.27 M Q 
yr" 1 for the SINFONI subsample of galaxies (their Ta- 
ble 3). In this scenario an initial, short, but massive, 
infall of gas creates a turbulent disk in which significant 
star formation occurs. After the resulting supernovae 
explosions clear the intercloud medium, there remains a 
massive, geometrically thick, collisional disk of dense gas 
clouds. Based on the observed properties, the Seyferts in 
the SINFONI subsample are thought to be in this second 
phase, which has a lifetime of at least 100 Myr. In this 
model, it is assumed that the physical properties of the 
gas on these scales are determined chiefly by the mass 
accretion rate from scales of ~100 pc. The disk height of 
the collisional disk is then determined by the cloud sizes, 
the transparency of the gas, and the cloud volume filling- 
factor, each of which can be estimated from the observed 
H2 kinematics and an assumed cloud mass of 10 Mq and 
a sound speed of 1.5 km s" 1 within the cloud. While 
the mass inflow rate is significant during this phase, it 
is slowly decreasing, and once it drops sufficiently the 



collisional disk transitions to become geometrically thin. 

In conclusion, on a scale of tens of parsecs there is no 
evidence that outflows or warps contribute significantly 
to the H2 1-0 S(l) kinematics of the sample galaxies, 
which would tend to exclude models invoking such mech- 
anisms to obscure the AGN at these radii. Although 
radiation pressure from AGN, stellar winds, and super- 
novae, are all likely to contribute to the H2 disk height 
at some level, none are capable of producing the disk 
heights at the radii considered here. At the peak of 
the star formation in the nuclear region, it is possible 
that radiation pressure from the stars can thicken the 
gas disk to the heights implied by the kinematics. How- 
ever, without knowing the evolution of the star formation 
with time it is not possible to say if the current star for- 
mation is capable of maintaining the vertical structure. 
Based on a model in which the gas properties on scales 
of tens of parsecs are dictated by the rate of inflow, the 
disk heights in the observed Seyferts may be maintained 
by energy input from the accretion of material into the 
nuclear region. Only with improved constraints on the 
physical structure of the interstellar medium on these 
scales from future higher spatial resolution observations, 
and increased sample sizes, will it be possible to quan- 
tify the relative contributions of the various mechanisms 
discussed to maintain the disks heights implied in the 
Seyfert galaxies observed. 

6. RELEVANCE TO AGN 

6.1. Are We Seeing an Obscuring Torus? 

In half of the sample galaxies (four of the eight), the 
extinction of the nuclear star formation derived in § 13.4.21 
implies gas column densities greater than the 10 22 cm" 2 
required to obscure the optical and UV emission from 
the AGN and its associated BLR (jRisaliti et alj I1999L 
iTreister et alj |2004). Of note is that the Seyfert 2 galaxy 
Circinus (and NGC 1068, which has been excluded from 
the analysis of H2; see §[3]) falls within this group of ob- 
scured AGN. In another two galaxies, the derived column 
densities are Njj ~ 10 22 cm" 2 , but considering the esti- 
mated errors it is uncertain if the densities are enough 
to obscure AGN emission. The column densities of the 
remaining two galaxies are below the threshold for ob- 
scuring the AGN. The additional clumpy component of 
dense gas implied by the high gas masses derived from the 
H2 kinematics (§ I3.4[) is, as discussed in § 13.51 unlikely to 
significantly contribute to obscuration of the AGN since 
only a few percent of the lines of sight will encounter a 
gas clump. We therefore conclude that in approximately 
half of the Seyfert galaxies observed the column density 
along lines of sight that pass through the molecular gas 
detected on scales of tens of parsecs will be great enough 
to obscure the central AGN. 

As discussed in § 13.31 and § the high velocity dis- 
persion and radial profile of the molecular gas on the 
scales observed suggests that the molecular gas disk 
in Seyfert galaxies is vertically thick due to bulk mo- 
tion of the gas, and thus some fraction of the lines of 
sight will pass through this disk. Observations of the 
space density of Seyfert Is (unobscured) versus Seyfert 
2s (obscured) suggest that the half-opening angle of the 
torus pre dicted by popular AGN unification schemes is 
39 °-48 (lOsterbrock fc MaTtdl fl99l iRush et all fl996l 
ISchmitt et al.ll2001llHao et alj|2005j) 7 implying that z Q /r 
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~ 0.9-1.3. The disk heights derived in £13.31 are consistent 
with this prediction with an estimated average z /r of 1.3 
±0.2. This fraction is unchanged if only those galaxies 
for which Njj > 10 22 cm~ 2 are considered. Therefore, if 
the molecular disks observed on scales of tens of parsecs 
are capable of obscuring AGN emission (as half of the 
galaxies observed suggest), then the derived disk heights 
imply the fraction of Seyferts observed through the thick 
obscuring disk is consistent with the observed space den- 
sities of Seyfert Is and 2s. 

Since the observed properties of the nuclear molecu- 
lar gas suggest a rotating, geometrically thick, clumpy, 
gas disk, we conclude that in Seyfert galaxies the gas 
on scales of tens of parsecs is likely to be associated 
with the smaller scale torus invoked in AGN unifica- 
tion models. This clumpy disklike structure is con- 
sistent wi th the radii predicted by models of clumpy 
tori (e.g. ICamerqn et al.lll993l . ISchartmann et al.ll2005l. 



Honig et all 120061 iFritz et all [20061. ISchartmann et al l 



20081 ISchartmann et al.ll2009f ). Given the estimated col- 



umn densities, it is also possible that the gas on these 
scales plays a role in obscuring the AGN. While there is 
evidence that obscuring material exists at smaller scales 
than the ~20 pc scales probed here (e.g. |Greenhill et al.l 
120031 Uaffe et al.l [200l iTristram et al.N2007t ). we would 



argue that the molecular gas measured reveals the larger 
scale extended structure of the nuclear obscuring mate- 
rial. Although it is most likely that the smaller scale 
torus component is largely responsible for obscuring the 
AGN, and other functions such as collimation of ioniza- 
tion cones, the contribution of the extended structure 
to AGN obscuration can be significant. Furthermore, 
clumpy torus models predict similar additional larger 
scale structures. For example, ISchartmann et all (2009) 
find that a clumpy torus consisting of a dense compact 
(~5 pc) component embedded in a more extended and 
diffuse component (~50 pc) is required to match the ob- 
served spectral energy distributions of Seyfert galaxies. 
In addition, the alignment of the clumpy disklike struc- 
ture with the large-scale orientation of the host galaxy 
disk is consistent with evidence that Seyfert galaxies are 
likely to fuel their AGN via secular processes rather than 
by the major mergers that are thought to trigger the fuel- 
ing of higher luminosity quasars (e.g. iDe Robertis et all 
[19981 iGarria-Burillo et alll2005h . 

6.2. Starburst-AGN connection 

As shown by a detailed analysis of the SINFONI 
subsample of Seyfert galaxies in D07, as well as 
by other studies of mostly Seyfert 2 galaxies (e.g . 
Gonzalez Delgado et al1l2001llCid Fernandes et~aTll2004l . 
Riffe l et all 12007ft . recent nuclear star formation is asso- 
ciated with the AGN activity. Since the star formation 
measured on scales similar to the H2 distribution is likely 
to have been Eddington limited (D07), it is inevitably 
short-lived. This implies that the star formation may 
be episodic in nature with multiple short bursts, the fre- 
quency of which would be dependent on the variability 
of the inflow of gas to the nuclear region (central ^100 
pc). With the typical velocity, velocity dispersion, and 
M gas values observed (e.g. 75 km s -1 , 50 km s" 1 , and 
10 7 M©, r espectively), an inflow rate of ~10 M Q yr _1 is 
expected ([Forster Schreiber et al.ll2006f ). This is reason- 
able for a short time period, but is certainly not likely to 



be sustainable. 

The correlation of the gas velocity dispersion with star 
formation that is shown in Fig. [5]suggests that the height 
of the molecular gas disk is related to the state of the 
nuclear star formation. As discussed in § (5[ the disk 
height may be maintained by stellar r adiation pressure 
(|Thompson. Quataert. &; Murravil2005f) . It is also possi- 
ble that the nuclear star formation and disk height are 
both dictated by the inflow rate i nto the nuclear region 
(|Vollmer. Beckert. fc Da vics 2008]) and that stellar radi- 
ation pressure is a secondary effect. This later scenario is 
consistent with the clumpy nature of the gas (§ 13. 5ft . since 
on these scales the gas is predicted to be in a geometri- 
cally thick, collisional disk that is dominated by compact 
clouds rather than a diffuse medium, and is thus trans- 
parent. Only immediately after a short, massive inflow 
of material into the central ^100 pc would the struc- 
ture be opaque, and this p hase lasts for only ~10 Myr. 
iGarcia-Burillo et"aH (|2008l ) suggest a possible triggering 
mechanism for this inflow event via the dynamical de- 
coupling of a nuclear bar with respect to an outer bar, 
which results in significant gas inflow. In addition, they 
find in numerical simulations that gravitational torques 
combined with viscosity can result in recurrent episodes 
of inflow. 

If the disk height is related to the state of the nuclear 
star formation and/or the inflow rate to the central ~100 
pc, we may expect that the thick molecular gas disks de- 
tected in the observed Seyferts is a dynamic structure. 
Each burst of gas inflow or star formation could increase 
the a of the gas, resulting in a thickening of the disk. 
In contrast, when the inflow rate decreases and nuclear 
star formation is in a low intensity state, the a of the 
gas would quickly decrease, allowing the gas to settle 
into a thin disk. In cases where the gas disk on these 
scales is capable of obscuring the AGN, as is the case 
in at least half of the observed Seyferts, a high state 
of inflow or star formation would result in an obscuring 
torus-like structure on scales of tens of parsecs and in 
a low state no obscuring medium on these scales would 
be present. In this case, the fraction of Seyfert Is ver- 
sus 2s would depend not only on orientation but also on 
the nuclear inflow rate and the level of star formation. 
This would imply that on average type 2 Seyferts have 
a higher rate of inflow into the nuclear region and/or 
a higher SFR per unit area out to radii of tens of par- 
secs. Unfortunately, it is not possible to rigorously test 
either of these predictions with the data currently avail- 
able. The star formation properties on scales of tens of 
parsecs have been investigated in only a small number 
of Seyferts (e.g. D07), and the inflow rate into the nu- 
clear region has been determined in even fewer AGNs, 
the majority o f which are at most low lum inosity Seyfert 
galaxies (e.g. iGarria-Burillo et all 12001 . HST imag- 
ing studies have however, found that the central kpc of 
Seyfert 2 galaxies exhibit a more disturbed morphology 
(more twisting of the isophotes) than is seen in Seyfert 1 
galaxies, wh ich can be interpreted as a sign of a greater 
inflow rate ([Hunt fc Malkanl I2004D . It is important to 
note that obscuration of the AGN is also likely to occur 
on scales smaller than those probed here. Therefore even 
if a state of low inflow and/or star formation leads to a 
thin disk on scales of tens of parsecs, and thus less prob- 
ability of obscuration of the AGN, it is still possible for 
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Seyfert 2-like properties to be observed. 

7. CONCLUSIONS 

With a mean spatial resolution of ^20 pc (ranging from 
4 to 36 pc) we have mapped the 2-D distribution and 
kinematics of H2 in the central ~100 pc of a sample of 
local AGN, most of which are Seyfert 1 galaxies (i.e. have 
a direct view of the nuclear region). From this sample we 
conclude that, with few exceptions (notably NGC 1068), 
the molecular gas is in a rotating disklike distribution 
that is geometrically thick due to bulk random motion. 
The HWHM of this disk is ~30 pc and the kinematics, 
as well as distribution, of the molecular gas is consistent 
with that of the stars in the nuclear region, suggesting 
that the stars and gas are spatially mixed on these scales. 

The disk height implied by the velocity dispersion of 
the gas is on the order of the disk radius [z ~ r). The gas 
column density is determined using the typical gas mass 
fraction of f g ~ 10% suggested by four independent esti- 
mates. This value implies an average Njj= (4.9±3.3) x 
10 23 cm -2 for the sample of galaxies, with values ranging 
from 1.3 to 10.0 x 10 23 cm~ 2 . However, the estimated 
stellar extinction within a similar region implies an aver- 
age gas column density of only Njj = (2.4±1.6) x 10 22 
cm~ 2 . The molecular gas on these scales is therefore in 
a clumpy distribution rather than uniformly dispersed. 
Although in some Seyfert galaxies the column densities 
derived by the stellar extinction are great enough to ob- 
scure an AGN, in approximately half of the galaxies the 
gas on these scales is transparent despite the ~10 7 M© 
of gas present within the central ^60 pc. We conclude 
that the molecular gas on the ~20 pc scales observed in 
the sample Seyfert galaxies is in a generally rotating, ge- 
ometrically thick clumpy disk that is associated with the 
smaller scale obscuring torus invoked in popular AGN 
unification models, but is not in all cases obscuring the 
AGN. 

Given the time-averaged SFRs and the estimated peak 
SFRs for the SINFONI subsample of galaxies, the verti- 
cal height of the gas disk can be produced through ra- 
diation pressure from the nuclear star formation occur- 
ring within the nuclear region. However, it is uncertain 
if the current SFRs are able to maintain the full disk 
height implied by the observed velocity dispersion. Ra- 
diation pressure from the AGN, stellar winds, and su- 
pernovae are also likely to play a role, albeit at a less 
significant level. There is no evidence that outflows or 



disk warps are a general mechanism for obscuring the 
AGN in these Seyfert galaxies, as there is no indica- 
tion of either in the kinematics measured down to the 
smallest scales. Alternatively, the disk heights may be 
maintained via converted gravitational energy as a re- 
sul t of gas accretion into the nu c lear r egion, as suggested 
bv IVollmer. Beckert. fc Daviesl (|2008f) . The kinematics 
of the nuclear gas disks within the central ~100 pc are 
consistent with the larger scale orientations of the galaxy 
disks in the sample galaxies, and, for those galaxies for 
which single-slit ISAAC spectra were obtained, the small- 
scale kinematics are also consistent with the rotation of 
the galaxies on kpc scales. 

Since the vertical structure of the gas disk may be de- 
pendent on the inflow rate and/or intensity of the nuclear 
star formation, we propose that the thick disks observed 
may be a dynamical structure. Significant nuclear gas in- 
flow or star formation will form a thick disk, while with 
low levels of inflow and/or star formation the gas will 
quickly cool into a thin disk. This scenario is supported 
by the correlation of the time-averaged SFR per unit 
area with the H2 velocity dispersion, and by the consis- 
tency of the nuclear star formation and gas densities with 
"Schmidt-Kennicutt" relations observed in other galax- 
ies. For those galaxies in which the gas on scales of tens 
of parsecs plays a role in obscuring the AGN (approxi- 
mately half of the observed galaxies), then this scenario 
suggests the nuclear star formation and/or gas inflow 
rate will in part determine the classification of a Seyfert 
galaxy as type 1 or 2. 
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APPENDIX 

A. DESCRIPTION OF HOT MOLECULAR GAS IN INDIVIDUAL OBJECTS 

A.l. NGC 1097 

NGC 1097 has an H2 flux distribution that is centrally concentrated with an elongation in a roughly north-south 
direction and a peak flux coincident with both the JT-band continuum and nonstellar emission. The distribution has 
an HWHM of 28.5 pc, which is resolved based on a spatial resolution determined from the nonstellar emission of 0".25, 
or 22 pc at the distance of NGC 1097. A Sersic fit to the azimuthal average of the flux distribution gives a disklike 
fit with 7i=1.5. The smoothness of the H2 velocity field indicates there is no significant warp in the gas disk down to 
the smallest scales measured. The best-fit kinematic parameters to the gas velocity field are consistent with that of 
the stars, with a best fit PA of -52 and an inclination of 42 given by the gas and PA of -49 and an inclination of 
43 by the stars. The kinemetry fit to the gas velocity field reveal deviations from co-planar disk rotation of the order 
of 30 km s _1 (greater than the ~10 km s _1 error of the velocity measurements) and outline a spiral structure. An 
investigation of the detailed kinematics in NGC 1097 will be presented in a future paper. 

As reported by D07, the central stellar velocity dispersion decreases significantly, displaying what has become known 
as a "sigma-drop" . Within the central 1" the dispersion drops to 100 km s _1 compared to a relatively constant value 
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of 150 km s -1 beyond this radius. Based on this sigma-drop, and evidence of recent star formation, the presence 
of a dynamically cold nuclear disk is suggested. The kinematics of the gas supports this in that the central a and 
rotational velocity change within the central 1" to more closely match that measured for the stars. The gas a increases 
from an average of 45 km s -1 outside of r = 0".5 to an average of 63 km s _1 and a maximum of 74 km s _1 at the 
nucleus (Fig [5] and Table [4j> . The offset of about 10 km s _1 higher a measured by ISAAC compared to that measured 
with SINFONI is due to blurring of the steep velocity gradient (~300 km s _1 in the central 2") in the lower spatial 
resolution seeing-limited ISAAC data. 

The dynamical mass derived from the gas kinematics is Md yn — 9.5 x 10 7 M at 30 pc, and is similar to the 
HWHM radius of 28.5 pc. This mass is lower by a factor of 1.4 than the dynamical mass determined from the stellar 
kinematics (D07), but is consistent when the velocity and a errors of ±10 km s _1 are considered. The kinematics 
within this region is there fore dominated by the BH mass, which is estimated to be within the range of 4-12 x 10 7 M© 
(jLewis fc Eracleousll2006f ). Although the H2 velocity dispersion is not as significant in comparison to the rotational 
velocity as is seen in other galaxies in the sample {V ro t/<J reaches unity at a radius of 24 pc, only slightly greater that 
the spatial resolution), the disk height at a radius of 30 pc is still estimated to be comparable to the radial scale. The 
gas column density based on the Md yn and the adopted f g — 10% is above 10 23 cm -2 out to over 150 pc, with Njj = 
4.2 x 10 23 cm -2 at 30 pc. Based on the cold molecular gas measured in CO (|Hsieh et al.ll2008| ). the gas mass fraction 
may be as low as 1.3% (Tabled, which lowers the column density to Njj = 5.5 x 10 22 cm -2 at 30 pc. The column 
density implied by the extinction of the stellar continuum based on the mixed model is 1.7 x 10 22 cm which is just 
above the 10 22 cm -2 limit of obscuring the AGN. 

A.2. NGC 3227 

The flux distribution of the hot molecular gas in NGC 3227 is relatively complex compared to the other galaxies in the 
sample. The strongest H2 TO S(l) emission is not found at the location of the nucleus, but is instead distributed along 
a PA o f -45° with peaks of emission 0".24 northwest and 0".35 and 0".81 southeast of the nucleus. ISchinnerer et al.l 
(2000a) also report that strongest emission from the cold molecular gas, measured by CO on a larger scale of 0".6, is 
not coincident with the nucleus. The azimuthal average of this complex emission structure peaks at a radius of about 
0".13. The azimuthal average of the OSIRIS data differs from that measured with SINFONI because of the FOV 
coverage (Fig. [TTj) . The OSIRIS data only covers east of the nucleus, and therefore does not include the northwest 
emission knot, which largely is responsible for the off-center emission peak of the full FOV azimuthal average measured 
with SINFONI. For this reason, the properties derived from the SINFONI data are taken to be more representative of 
the full nuclear region in NGC 3227. 

With an azimuthal average it is of course not possible to represent the complex emission structure seen in NGC 
3227, and a Sersic fit to the averaged distribution at the radii of the nonsymmetric emission is thus of limited value. 
We therefore fit a Sersic function to the azimuthally averaged flux distribution beyond the radii effected by the knots 
of emission, and the HWHM of the observed distribution can then be estimated from an extrapolation of this fit. 
For the SINFNOI data a fit is preformed to the fall in emission beyond the peak at 0".13 (10 pc), and a best fit is 
found with n=1.76. The OSIRIS data are fit outside a radius of 0".39 (30 pc), which is beyond the majority of the 
bright emission knots and most likely represents the more extended emission. The best fit in this case is found with 
a n=1.16 Sersic function. The HWHMs of the extrapolated Sersic fits are 29.3 pc and 34.9 pc for the SINFONI and 
OSIRIS data, respectively. For comparison, the HWHM of the azimuthally averaged observed distribution measured 
with OSIRIS is 60.9 pc with respect to the nucleus and 28.1 pc with respect to the peak of the distribution. These 
HWHM estimates, both the observed values and those determined from the Sersic fits, are well resolved based on the 
spatial resolution given by the broad Br7 emission, which for the SINFONI data is 0".085 (7 pc) and for the OSIRIS 
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Fig. 11. — Comparison of the azimuthally averaged (a) flux distribution and (b) velocity (triangles) and velocity dispersion ("x" symbols) 
in NGC 3227 measured with SINFONI and OSIRIS. As indicated by the legend, the solid lines are for SINFONI measured data and 
dash-dotted lines are those measurements made with OSIRIS. The flux distribution has been normalized to the peak of the extrapolated 
Sersic fits, which for the S INFO NI (solid curve) and OSIRIS (dashed curve) data were fit beyond a radius of 10 pc (0".13) and 31 pc 
(0".39), respectively (see § IA.2l for details). 
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data is 0".07 (6 pc). 

Despite the complexity of the H2 flux distribution, the kinematics of the molecular gas is fairly regular. The best fit 
inclination and PAs of the gas disk determined from the SINFONI and OSIRIS data are consistent, with i — 50°-60° 
and PA = -30 to -45 . These disk parameters are in agreeme nt with the best fit to the st ellar kinematics o f i = 55 °-60 ° 
and P A = -45 , as well as with the large-scale galactic disk (jde Vaucouleurs et al.lll99ll ) . As discussed in iDavies eT ID 
I2006L the majority of the gas is in uniform rotation, with only the weak emission along the minor axis showing signs 
of nonrotational motions, possibly an outflow. The kinemetry fit shows deviations from co-planar disk rotation along 
the minor axis on the order of ~30 km s" 1 (compared to the ~10 km s _1 error of the velocity measurements). The 
velocity dispersion is rather high, with an average value of ~86 km s _1 within r = 0".5 (94 km s _1 measured with 
SINFONI and 77 km s _1 measured with OSIRIS). The velocity field, as well as the velocity dispersion, is comparable 
to that measured for the stars in the nuclear region, with a maximum difference of 25 km s _1 along the minor axis. 
A comparison of the azimuthally averaged velocity and velocity dispersion as measured with SINFONI and OSIRIS is 
shown in Fig. 1111 

The gas kinematics indicate that M& vn = 20 x 10 7 M at a radius of 30 pc. Modeling of the H2 1-0 S(l) kinematics 
suggests that NGC 3227 has a BH of 2.0± x 10 7 M (Hicks fc Malkanll2008h . and modeling of the stellar kinematics 
gives a mass estimate of 0.7-2.0 x 10 7 M Q (|Davies et al.ll2006[) . therefore the BH is not a significant fraction of the 
dynamical mass measured. The H2 velocity dispersion is greater than the inclination corrected rotational speed out 
beyond 30 pc, suggesting that the gas is in a thick disk ~45 pc high at 30 pc from the nucleus. Although measurements 
of the cold molecular gas detected with CO (2-1) show that V ro t/cr greater than one outside of 1", even the cold gas 
has a significant velocity dispersion of 60 km s _1 inside of this radius (Table [5J ISchinnerer et al.ll2000af ). Assuming 
f g = 10% (which is consistent with estimates derived from CO observations; Table EJ, and that this gas is uniformly 
distributed, the average gas column density inside of 30 pc is estimated to be ~ 9 x 10 23 cm -2 . Based on the 
stellar extinction a mixed model suggests that the column density is 1.9 x 10 22 cm' 2 , which is well over an order of 
magnitude lower than that derived from the dynamical mass, suggesting that the majority of the nuclear gas is in a 
clumpy distribution. 

A more detail e d ana lysis of the stellar and gas kinematics in NGC 3227 is presented in IDavies et alJ (|2006l ) and 
iHicks fc Malkanl d2008h . 

A.3. NGC 3783 

The H2 1-0 S(l) flux distribution in NGC 3783 is centrally concentrated and the strongest emission is consistent 
with the nucleus location. The observations of NGC 3783 have a spatial resolution (based on both the nonstellar 
emission and the broad Br7 emission) of 0".18, which is equivalent to 36 pc at the galaxy distance. A Sersic fit to the 
azimuthally averaged light distribution gives a best fit with n = 1.85 and the azimuthally averaged light profile has an 
observed HWHM of 20.2 pc. 

The H2 kinematics can reliably be measured out to only about 1" (equivalent to 200 pc). The velocity field is 
relatively regular, with a best-fit PA of -15 °, and inclination of 35 °. Over the central 0".5 less than 4% of the if-band 
continuum is stellar (D07), which prevents a reliable measurement of the stellar kinematics, and thus no comparison 
with the gas kinematics is possible. A comparison to the if-band photometry (Mulchaey et al. 1997) indicates that 
the gas kinematics is consistent with the large scale isophotes, suggesting that there is no warp in the gas disk. 

At a radius of 30 pc, a dynamical mass of Md yn — 2.8 x 10 7 M Q is derived from the gas kinematics, and at 
the observed HWHM of 20.2 pc M dyn = 1.6 x 10 7 M Q . Based on the mass estimated from reverberation mapping 
(|Peterson et al.l |2004[ ) the central BH has a mass comparable to the dynamical mass within a radius of 30 pc and 
therefore it likely dominates the kinematics. The disk height is estimated to be similar to the radial scale, with a 
vertical scale of 40 pc at a 30 pc radius. The gas column density in NGC 3783 (assuming f g — 10%) is a bit lower than 
estimated for the other galaxies in the sample, but is still above 10 23 cm -2 out to a radius of ^70 pc and is 1.3 x 10 23 
cm -2 at 30 pc. The stellar extinction implies Njj = 5.3 x 10 22 cm -2 , which is still high enough to obscure the AGN. 

A.4. NGC 4051 

In NGC 4051, the H2 flux distribution is symmetric and centrally concentrated at a location coincident with the 
nonstellar emission. The observed azimuthally averaged HWHM is 7.6 pc, which is just barely resolved based on the 
spatial resolution determined from the nonstellar emission of 0".12, or 6 pc at the distance of NGC 4015. The best fit 
to the azimuthally averaged flux distribution is a n=1.8 Sersic function. 

The measured velocity field is consistent with uniform disk rotation without any indication of a warp down to the 
smallest scales measured, but the limited FOV and relatively low signal-to-noise ratio of these data do not constrain 
the kinematics as well as in the other galaxies. The best fit to the H2 kinematics suggests a disk with i = 50° and 
PA = 5 °. At a radius of 30 pc the V ro t/cr ratio is 0.83 and a disk height of z a = 41 pc is derived from the kinematics. 
The dynamical mass within this rad ius is 5.0 x 10 7 M and, based on the reverberation mapped Mbh estimate of 
0.2 x 10 7 M (jPeterson et al.l 12004). the Mbh does not dominate the kinematics on these scales. The derived gas 
column density within r = 30 pc (assuming f g = 10%) is then 2.2 x 10 23 cm -2 and remains above 10 23 cm -2 within 
the measured FOV. Unlike all other galaxies in the sample, there is very little extinction of the nuclear stars with a 
mixed extinction model predicting at most Ay = 1 mag. 
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A.5. NGC 4151 

The H2 flux distribution in the nuclear region of NGC 4151 is relatively complex due to knots of emission along 
a PA of roughly 90°. As noted bvlHicks fe Malkan ] (I2008D . several of these knots are coincident with radio emission 
associated with a nuclear jet (e.g. iMundell et al.l l2003). Underlying these knots of emission is a component of extended 
H2 emission that is centrally peaked at the location of the nonstellar emission. We do not attempt to examine the 
bright knots of emission further in this paper, and exclude them from the azimuthal average by only considering 
flux below 25% of the peak. Furthermore, H2 is not well detected in the central 0".2, and this region is binned to 
obtain a reliable measure of the central flux and kinematics. Unfortunately, a treatment of the bright emission knots 
such as that carried out for NGC 3227 is not possible because the knots extent throughout the measured FOV. The 
azimuthally averaged flux distribution of the extended H2 emission has an observed HWHM of 28.1 pc and is best fit 
with ann= 1.3 Sersic function. If instead all of the emission is considered in the azimuthal average then there are 
two additional peaks at ~0".4 and 0".7 due to knots of emission on either side of the nucleus at these radii. In Fig. 0] 
this radial profile (smaller diamonds with dashed error bars) can be compared with that of the diffuse emission (solid 
diamonds and error bars). Based on the nonstellar emission, the spatial resolution of these NGC 4151 data is 0".ll, 
or 7.5 pc at the distance of NGC 4151. 

The underlying H2 emission has a velocity field suggestive of uniform disk rotation with a best fit of PA = 10 and i 
= 45 °, which is consist ent with the kinematics of th e stars in the nuclear region and roughly consistent with the large 
scale galactic rotation (jde Vaucouleurs et al.l[l99ll ). Coincident with the knots of bright emission mentioned above 
are patches of velocity inconsistent with co-planar circular rotation, with residuals on the order of ^30 km s . This 
deviation from rotation is not surprising if indeed the knots of emission are associated with the radio jet. The velocity 
dispersion throughout the measured FOV is ~55 km s , and is greater than the inclination corrected rotational 
velocity out to a radius of about 20 pc, beyond which V ro t/<J rises steadily above unity. The kinematics at this radius 
implies a disk height of z a = 33 pc and a dynamical mass of 14.0 x 10 7 Mfq, suggesting that the central BH does no t 
dominate the kinematics at this radius based on the estimated BH mass ( Hic ks fc Malkanll2008l iPeterson et"aTll2004h . 
The gas column density remains above 10 23 cm~ 2 out to the edge of the FOV and is estimated to be 6.2 x 10 23 cm~ 2 
at a radius of 30 pc, assuming f g = 10%. A column density of an order of magnitude less is predicted by fitting a 
mixed extinction model to the stellar continuum. The extinction suggests that Nh = 2.8 x 10 22 cm -2 , which is still 
high enough to obscure the AGN. 

A.6. NGC 6814 

The H2 emission detected in NGC 6814 is nearly symmetric and peaks at the AGN location of the nonstellar emission. 
With a spatial resolution of 0".17 (18 pc), the observed HWHM of the azimuthally averaged distribution is measured 
to be 34.8 pc. Fitting a Sersic function indicates a disklike distribution with n = 1.04. 

Like the other galaxies in the sample, NGC 6814 exhibits regular rotation indicative of uniform disk rotation. The 
best fit to the gas kinematics suggests a disk with i = -40 and PA = 55 °, which is consistent with results from fitting 
the stellar kinematics. The velocity dispersion is greater than the inclination corrected rotational velocity out to the 
point where a reliable fit can no longer be obtained at about 0".5. At a radius of 30 pc, the disk height based on the 
kinematics is estimated to be comparable to the radial distance with z Q — 46 pc. The dyn amical mass a t this radius is 
4.2 x 10 7 M©, a quarter of which is estimated to be mass attributable to the central BH (Wandcl 2002). The derived 
gas column density is greater than 10 23 cm -2 to a radius of 40 pc, with Njj= 1.9 x 10 23 cm" 2 at 30 pc assuming f g = 
10%. The stellar extinction in this galaxy is relatively low, with an implied column density of only Njj — 5.7 x 10 21 
cm~ 2 . This estimate is well over an order of magnitude lower than derived from the dynamical mass, suggesting that 
the majority of the gas in this galaxy, on scales of tens of parsecs, is in a clumpy distribution. 

A.7. NGC 7469 

The H2 1-0 S(l) flux distribution in NGC 7469 is very symmetric with a peak coincident with both the nonstellar and 
stellar emission. The spatial resolution of the SINFONI data is estimated from the nonstellar continuum to be 0".15, 
which is 48 pc at the distance of NGC 7469. The azimuthal average of the SINFONI observed distribution has an 
HWHM of 77.4 pc, and a fit to the azimuthal average finds a Sersic n = 1.3 distribution. The OSIRIS data has a higher 
spatial resolution of 0".ll, equivalent to 35 pc, and for this reason the properties measured from the OSIRIS have data 
are taken to be more accurate, although, with the exception of the observed flux HWHM, the measured and derived 
properties from both data sets are consistent. The observed HWHM of the azimuthally averaged flux distribution 
measured with OSIRIS is 30.8 pc and a best fit to the flux profile is found with a Sersic n = 1.98 distribution. A 
comparison of the azimuthally averaged flux distribution measured with OSIRIS and SINFONI, as well as that of the 
velocity and velocity dispersion, is shown in Fig. 1121 

The H2 velocity field is suggestive of uniform rotation, with no indication of a significant warp in the gas disk down 
to the smallest scales measured. The gas kinematics is best fit with i — 50 and PA = -40 in the SINFONI data, in 
comparison to i = 30 and PA = -48 in the higher spatial resolution OSIRIS data. These results just barely agree 
to within the expected ±10° errors of the measure ments, and both agre e with the fit to the nuclear stellar kinematics 
and with the larger scale gas disk measured in CO (jDavies et al.l l2004b). The velocity dispersion is about 60 km s _1 
in the center and rises as high as 89 km s _1 at r = 0".5. This dispersion in velocity is greater than the (inclination 
corrected) rotational velocity in the central 0".29 and V ro t/<J remains near unity out to at least 1". The kinematics 
implies a disk height of about 43 pc at a radius of 30 pc. The dynamical mass derived from the gas kinematics gives 
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Mdyn = 10.1 x 10 7 Mq within r = 30 pc, and assuming f g = 10%, the derived gas column density within 30 pc is 
4.2 x 10 23 cm~ 2 . Based on CO observations of the cold molecular gas (Tabled Davies eTaLll2004bn . the gas mass 
fraction is likely to have a significantly higher value of ~60%, suggesting an average gas column density of 2.4 x 10 24 
cm~ 2 out to a radius of 30 pc. A mixed extinction model also suggests column densities high enough to obscure the 
AGN, with Nh ~ 5 x 10 22 cm -2 . This column density is still an order of magnitude or more lower than implied 
by the dynamical mass, suggesting that a significant fraction of the nuclear gas is distributed in clumps rather than 
uniformly. 

A. 8. Circinus 

The H 2 1-0 S(l) emission in Circinus is symmetric with a peak coincident with that of the JT-band continuum. The 
HWHM of the azimuthally averaged light profile is 7.4 pc, observed with a spatial resolution of 0".22, or 4.2 pc at the 
galaxy distance. The light profile is best fit with a n = 1.1 Sersic function. The H2 velocity field suggests uniform 
rotation, with the exception of a region at the western edge of the FOV at 9 pc, which is most likely associated with the 
ionization cone. A gas disk is best fit to the kinematics with a PA of 30 and an inclination of 55 °. Wi thin the errors, 
these parameters agree with the large scale optical photometry (PA = 20° and inclination of 65°; Fre emen et al.l 
Il977f) . suggesting t hat there is no warp d own to at least a few parsecs (there is evidence of a warp on much smaller 
scales of 0.1-0.4 pc; iGreenhill et~aT1l2003h . 

The dynamical mass is estimated to be Md yn — 1.9x 10 7 M Q within 9 pc, and Md yn = 1.8x 10 7 M Q within the 
HWHM radius of 7.4 pc. Based on this Md yn and assuming a gas fraction of 10%, the gas column density within 9 
pc is 9.4 x 10 23 cm -2 , and increases to 1.3 x 10 24 cm -2 within the HWHM radius. The stellar continuum is heavily 
obscured in this galaxy, consistent with its classification as a Seyfert 2. A mixed extinction model implies column 
densities of Ng — 4.3 x 10 22 cm -2 , which is high enough to obscure the AGN. This estimate is still over an order of 
magnitude lower than that derived from the dynamical mass, suggesting a clumpy medium. The velocity dispersion is 
relatively constant throughout the FOV measured (a ~ 56±3 km s _1 for r<0".5; Table HJ) and the rotational velocity 
to a ratio is relatively low with V ro t/u only reaching 0.7 at the edge of the FOV. At this radius, the disk height is 
estimated to be 13 pc, and a height of 11 pc is estimated at the HWHM radius. 

A detailed analysis of the nuclear reg ion of Circinus, including the H 2 properties as measured with SINFONI, is 
presented in lMiiller Sanchez et alJ (|2006h . 

A.9. NGC 1068 

NGC 1068 has the most complex H 2 1-0 S(l) properties of any of the galaxies in the sample and is discussed in 
detail in iMiiller Sanchez et aT] (|2009ft . Within the central 1", there are two knots of H 2 emission, one coincident with 
the nucleus and another to the north, and there is no detectable H 2 emission in the southwest. An azimuthal average 
of this distribution gives an observed HWHM of 13.2 pc and a Sersic fit finds n — 2.09, although, like in NGC 3227, 
these values obviously do not reflect the two dimensional nature of the nuclear H 2 distribution. The nuclear stars 
are heavily obscured in this galaxy, and even after correcting for contamination from the nonstellar continuum, the 
extinction is so great that it is beyond the saturation limit of a mixed extinction model (which is at an optical depth 
of a few). The spatial resolution based on the nonstellar emission is 0".08, equivalent to 5 pc at the distance of NGC 
1068. 

Noncircular motions dominate the H 2 kinematics, and Miiflcr San chez et al.l (|2009l ) argue that streaming of the gas 
toward the nucleus is suggested. This should therefore be kept in mind when evaluating the azimuthally averaged 
rotational velocity and a shown in Fig. [3l Because the H 2 properties in NGC 1068 so greatly differ from the disk 
rotation seen in the stars, unlike the consistency seen in the other sample galaxies, this galaxy is not included in the 
analysis of general H 2 properties in AGNs. This galaxy therefore represents an exception to the general conclusions 
presented. 




Fig. 12. — Comparison of the azimuthally averaged (a) flux distribution and (b) velocity (triangles) and velocity dispersion ("x" symbols) 
in NGC 7469 measured with SINFONI and OSIRIS. As indicated by the legend, the solid lines are for SINFONI measured data and 
dash-dotted lines are those measurements made with OSIRIS. 
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